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ABSTRACT 


An  evaluation  of  the  performance  of  the  coherent  sidelobe 
cancellation  (CSLC)  techniques  in  suppressing  nighttime  co-channel 
communications  interference  and  ionospheric  clutter  in  HF  surface  wave 
radars  operating  in  coastal  regions  was  carried  out.  In  particular,  the 
optimum  order  of  execution  among  the  three  signal  processing  functions: 
(i)  Doppler  processing,  (ii)  digital  beamforming  and  (iii)  interference 
cancellation  was  of  interest,  in  terms  of  signal-to-interference-plus-noise 
ratio  (SINR)  improvement  and  angular  accuracy.  The  optimum  order  of 
execution  between  the  above  mentioned  operations  depends  on  the 
characteristics  of  the  interference.  In  terms  of  ionospheric  clutter,  there  is 
insufficient  information  to  arrive  at  definitive  conclusions.  A  preliminary 
examination  of  the  temporal,  spatial  and  spectral  characteristics  of 
ionospheric  clutter  was  carried  out  to  aid  in  the  development  of  mitigating 
algorithms. 


RESUME 


On  a  evalue  la  performance  de  techniques  d'annulation  coherente 
des  lobes  lateraux  (CSLC)  pour  la  suppression  des  interferences  nocturnes 
des  communications  radio  dans  le  meme  canal  et  du  clutter  ionospherique 
des  radars  a  onde  de  surface  HF  exploites  dans  les  regions  coheres.  En 
particulier,  il  s'agissait  de  determiner  1'ordre  optimal  d'execution  de  trois 
fonctions  de  traitement  des  signaux,  soit  (i)  le  traitement  Doppler,  (ii)  la 
conformation  numerique  du  faisceau  et  (iii)  1'annulation  des  interferences, 
a  partir  de  l'amelioration  du  rapport  signal/interference  plus  bruit  (SINR) 
et  de  la  precision  angulaire.  L'ordre  optimal  d'execution  de  ces  fonctions 
depend  des  caracteristiques  des  interferences.  En  ce  qui  a  trait  au  clutter 
ionospherique,  il  n'existe  pas  suffisamment  d’information  pour  quWi  puisse 
tirer  des  conclusions  definitives.  Un  examen  preliminaire  des 
caracteristiques  temporelles,  spatiales  et  spectrales  du  clutter  ionospherique 
a  ete  effectue  dans  le  but  d'aider  a  etablir  des  algorithmes  d’attenuation. 
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EXECUTIVE  SUMMARY 


High  frequency  radars  rely  heavily  on  digital  signal  processing  to  improve  target-detection 
performance.  Central  to  the  signal  processing  suite  employed  in  HF  radars  are  Doppler  processing 
and  digital  beam-forming  operations.  In  a  practical  operating  environment,  however,  an  HF  radar 
must  be  able  to  function  in  the  presence  of  severe  interference  that  is  orders  of  magnitude  higher 
than  the  target  return  of  interest.  For  high  frequency  radars  operating  in  coastal  regions,  nighttime 
co-channel  communications  interference  and  direct  over-head  ionospheric  clutter  are  of  particular 
concern. 

Co-channel  communications  interference  is  the  dominant  source  of  interference  at 
nighttime.  Occasionally  it  is  present  during  day-light  hours.  Ionospheric  clutter  is  a  self-induced 
interference.  Its  effects  are  very  much  dependent  on  the  system  parameters  such  as  operating 
frequency,  antenna  pattern,  etc,  and  environmental  parameters  such  as  time-of-day,  sun-spot 
cycle,  etc. 

Because  both  of  these  interferences  possess  some  directional  characteristics,  they  can  be 
reduced  by  cancellation  techniques.  These  techniques  require  a  number  of  auxiliary  antennas,  the 
signals  from  which  provide  a  means  to  estimate  and  subtract  the  interference  from  the  main 
antenna  outputs.  These  are  collectively  referred  to  as  coherent  sidelobe  cancellation  (CSLC) 
techniques.  The  basic  structures  of  the  cancellation  schemes  used  against  ionospheric  clutter  and 
communications  interference  are  similar.  The  major  difference  lies  in  the  way  the  correlation 
properties  are  estimated.  The  process  of  interference  cancellation  modifies  the  gain  and  phase 
responses  of  the  array  elements.  A  natural  question  that  arises  is  whether  there  is  a  difference  in 
the  bearing  estimate  of  a  target  as  a  result  of  changing  the  order  of  execution  of  Doppler 
processing,  interference  cancellation  and  digital  beamforming.  The  work  carried  out  in  this  report 
is  aimed  at  obtaining  some  quantitative  answers  to  the  above  question. 

The  results  of  the  analysis  indicated  that  performing  CSLC  before  Doppler  processing 
yielded  the  best  result,  both  in  SINR  (signal-to-interference-plus-noise  ratio)  improvement  and 
in  angle-estimate  accuracy,  under  the  condition  of  a  single  dominant  interference  source. 
Performing  CSLC  after  Doppler  processing  yielded  the  best  result  when  the  requirement  is  to 
suppress  a  weaker  interference  signal  in  the  presence  of  a  more  dominant  interference  source 
such  as  sea  clutter  or  a  stronger  interference  signal  from  another  direction. 

Because  ionospheric  clutter  is  a  self-induced  interference  and  is  isolated  in  both  range  and 
Doppler,  certain  techniques  that  are  effective  with  communications  interference  in  extracting  the 
correlation  properties  among  the  main  and  auxiliary  array  antennas  are  ineffective  against 
ionospheric  clutter.  A  preliminary  investigation  of  the  temporal,  spatial  and  spectral  properties 
of  ionospheric  clutter  was  carried  out.  The  results  showed  that  the  ionospheric  clutter  does 
originate  from  a  finite  number  of  directions  and  there  appeared  to  be  short-term  stability  in  its 
correlation  properties.  Because  of  the  available  data  were  collected  from  a  linear  array,  it  was  not 
possible  to  systematically  evaluate  various  schemes.  Consequently,  it  can  only  be  conjectured  that 
time-domain  cancellation  is  probably  more  effective  than  Doppler-domain  cancellation,  although 
this  is  by  no  means  certain. 
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The  observed  characteristics  of  ionospheric  clutter  are  highly  dependent  on  the  operating 
radar  frequency  and  time-of-day.  The  strong  ionospheric  clutter  was  observed  to  migrate  in  range 
on  a  relatively  short  time  scale  (tens  of  minutes).  Because  the  range  extent  of  the  ionospheric 
clutter  is  relatively  narrow,  it  may  be  possible  to  mitigate  the  effects  of  ionospheric  clutter  by 
a  judicious  use  of  frequency  agility  and  advance  tracking  techniques. 

High  frequency  surface  radar  has  the  potential  to  be  an  important  wide  area  coastal 
surveillance  sensor  for  Canada.  The  results  of  this  work  contribute  to  the  continuing  effort  at 
DREO  to  improve  the  performance  of  this  type  of  radars.  Raytheon  Systems  Canada  Limited 
(RSCL)  has  been  an  industrial  partner  of  DND  in  the  research  and  development  of  HF  radars. 
This  work  was  partially  supported  by  RSCL  through  a  contract  with  DREO. 
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SOMMAIRE 


Les  radars  haute  frequence  reposent  largement  sur  le  traitement  numerique  des  signaux 
pour  ameliorer  la  detection  des  cibles.  Les  operations  de  traitement  Doppler  et  de  conformation 
numerique  du  faisceau  sont  essentielles  au  traitement  des  signaux  par  les  radars  HF.  Dans  un 
environnement  pratique,  un  radar  HF  doit  toutefois  etre  a  meme  de  fonctionner  en  presence 
d’ interferences  intense,  plusieurs  fois  superieures  a  l’echo  d’une  cible  d’interet.  Les 
interferences  nocturnes  des  communications  radios  dans  le  meme  canal  et  le  clutter 
ionospherique  revetentune  importance  particuliere  pour  les  radars  haute  frequence  exploites  dans 
les  regions  coheres. 

Le  principal  type  d’ interferences  durant  la  nuit  est  celui  des  communications  radios  dans 
le  meme  canal.  II  arrive  egalement  qu’il  se  produise  durant  le  jour.  Le  clutter  ionospherique 
est  un  type  de  d’ interferences  auto-induit.  Ses  effets  dependent  beaucoup  des  parametres  des 
systemes,  comme  la  frequence  de  fonctionnement  et  le  diagramme  d’antenne,  ainsi  que  des 
parametres  environnementaux,  comme  l’heure  et  le  cycle  des  taches  solaires. 

Possedant  tous  deux  certaines  caracteristiques  directionnelles,  ces  types  d’ interferences 
peuvent  etre  reduits  a  l’aide  de  techniques  d’annulation.  Ces  techniques  exigent  un  certain 
nombre  d’antennes  auxiliaires  dont  les  signaux  permettent  d’estimer  et  de  soustraire  les 
interferences  des  sorties  d’antenne  principales.  On  les  appelle  collectivement  techniques 
d’annulation  coherente  des  lobes  lateraux  (CSLC).  Les  structures  de  base  des  techniques 
d’annulation  utilisees  sont  semblables  pour  le  clutter  ionospherique  et  les  interferences  radios. 
La  principale  difference  reside  dans  la  fagon  dont  les  proprietes  de  correlation  sont  estimees. 
Le  processus  d’annulation  des  interferences  modifie  les  reponses  en  gain  et  en  phase  des 
elements  du  reseau.  La  question  se  pose  naturellement  de  savoir  s’il  existe  une  difference  de 
relevement  estime  d’une  cible  par  suite  d’un  changement  apporte  a  1’ordre  d’execution  du 
traitement  Doppler,  de  l’annulation  des  interferences  et  de  la  conformation  numerique  du 
faisceau.  Les  travaux  decrits  dans  ce  rapport  visent  a  obtenir  certaines  reponses  quantitatives 
a  la  question  ci-dessus. 

Les  resultats  de  l’analyse  ont  indique  que  l’execution  de  la  CSLC  avant  le  traitement 
Doppler  donnait  les  meilleurs  resultats,  tant  pour  ce  qui  est  de  1’ amelioration  du  rapport 
signal/interference  plus  bruit  (SINR)  que  de  la  precision  angulaire,  en  presence  d’une  seule 
source  d’ interference  dominante.  L’execution  de  la  CSLC  apres  le  traitement  Doppler  donnait 
les  meilleurs  resultats  lorsqu’il  fallait  supprimer  un  signal  d’ interference  plus  faible  en  presence 
d’une  source  d’ interference  dominante,  par  exemple  du  clutter  de  mer  ou  des  interferences  plus 
intenses  en  provenance  d’une  autre  direction. 

Comme  le  clutter  ionospherique  est  un  type  d’ interference  auto-induit,  isole  quant  a  la 
distance  et  a  1’effet  Doppler,  certaines  techniques  d’extraction  des  proprietes  de  correlation 
d’antennes  reseaux  principales  et  auxiliaires,  efficaces  pour  le  traitement  des  interferences  radios, 
s’averent  peu  performantes  contre  le  clutter  ionospherique.  On  a  effectue  une  etude  preliminaire 
des  proprietes  temporelles,  spatiales  et  spectrales  du  clutter  ionospherique.  Les  resultats  ont 
revele  que  le  clutter  ionospherique  ne  provient  pas  d’un  nombre  fini  de  directions  et  que  ses 
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configurations.  On  doit  done  se  limiter  a  supposer  que  1'annulation  dans  le  domaine  temporel 
est  probablement  plus  efficace  que  1'annulation  dans  le  domaine  Doppler,  bien  que  cela  ne  soit 
aucunement  confirme. 

Les  caracteristiques  observees  du  clutter  ionospherique  dependent  fortement  de  la 
frequence  de  fonctionnement  du  radar  et  de  l'heure.  On  a  observe  que  du  clutter  ionospherique 
intense  pouvait  se  deplacer  en  distance  sur  une  periode  de  temps  relativement  courte  (de  l'ordre 
des  dizaines  de  minutes).  Comme  les  limites  de  distance  du  clutter  ionospherique  sont 
relativement  etroites,  il  pourrait  etre  possible  d'attenuer  les  effets  du  clutter  ionospherique  en 
faisant  un  usage  judicieux  des  techniques  d'agilite  en  frequence  et  de  poursuite  anticipee. 

Le  radar  de  surface  haute  frequence  peut  devenir  un  important  outil  de  surveillance  de 
larges  regions  coheres  canadiennes.  Les  resultats  des  travaux  decrits  contribuent  aux  efforts 
permanents  du  CRDO  en  vue  d'ameliorer  les  performances  de  ce  type  de  radar.  Raytheon 
Systems  Canada  Limited  (RSCL)  est  un  des  partenaires  industriels  du  MDN  en  matiere  de 
recherche  et  developpement  de  radars  HF.  RSCL  contribue  partiellement  a  ces  travaux  dans  le 
cadre  d'un  marche  conclu  avec  le  CRDO. 


Chan,  Hing  C.,”Une  etude  sur  la  cancellation  d'interference  pour  les  radars  decametrique  a  onde 
de  surface,”  CRDO  TR2000-028,  Centre  de  Recherches  et  Developpement  Ottawa. 
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1. 


INTRODUCTION 


High  frequency  radars  rely  heavily  on  digital  signal  processing  to  improve  target-detection 
performance.  Central  to  the  signal  processing  suite  employed  in  HF  radars  are  Doppler  processing 
and  the  beam-forming  operations.  In  a  practical  operating  environment,  however,  an  HF  radar 
must  be  able  to  function  in  the  presence  of  severe  interference  that  is  orders  of  magnitude  higher 
than  the  target  of  interest.  To  minimize  the  detrimental  effects  of  interference,  cancellation 
techniques  are  used  to  exploit  the  directional  characteristics  of  the  interference  signals.  These 
techniques  require  a  number  of  auxiliary  antennas,  the  signals  from  which  provide  a  means  to 
estimate  and  subtract  the  interference  from  the  main  antenna  outputs.  This  subtraction  modifies 
the  gain  and  phase  responses  of  the  array  elements.  A  natural  question  that  arises  is  whether  there 
is  a  difference  in  the  bearing  estimate  of  a  target  as  a  result  of  changing  the  order  of  execution 
of  Doppler  processing,  interference  cancellation  and  digital  beamforming.  The  work  carried  out 
in  this  report  is  aimed  at  obtaining  some  quantitative  answers  to  the  above  question. 

For  high  frequency  radars  operating  in  coastal  regions,  there  are  four  major  sources  of 
interference,  namely,  sea  clutter,  ionospheric  clutter,  meteor  clutter  and  co-channel 
communications  interference.  For  sea  clutter,  there  have  been  many  attempts  at  its  cancellation. 
Most  of  the  techniques  that  have  been  reported  [1,2]  attempt  to  suppress  the  sea  clutter  in  the 
Doppler  domain.  At  the  present  time,  there  is  no  consensus  as  to  what  the  most  effective 
technique  is  regarding  sea-clutter  suppression.  Consequently,  sea  clutter  will  not  be  included  in 
this  study. 

Meteor  clutter  is  a  transient  phenomenon  in  that  the  duration  of  the  echoes  from 
meteoroid  trails  is  short  compared  with  the  coherent  integration  time  (CIT)  interval.  Because 
meteoroids  are  isolated  both  temporally  and  in  range,  it  is  very  difficult  to  obtain  reliable 
correlation  properties  of  the  array  samples  before  the  interference  vanishes.  A  reasonably 
effective  means  of  reducing  the  effects  of  meteor  clutter  is  impulse  excision. 

For  the  remaining  two  interference  sources,  there  are  effective  cancellation  techniques  that 
can  be  employed.  The  basic  structures  of  the  cancellation  schemes  against  ionospheric  clutter  and 
communications  interference  are  very  similar.  The  major  difference  lies  in  the  way  the  correlation 
properties  are  estimated. 

Ionospheric  clutter  is  a  self-generated  interference.  The  ionosphere  is  a  region  of  ionized 
gas  layers  (D,  E  and  F)  above  the  surface  of  the  earth.  The  echoes  from  these  ionized  layers 
comprise  both  a  specular  and  a  diffuse  component.  The  specular  component  is  confined  to  a 
small  number  of  range  and  Doppler  bins.  The  diffuse  component,  on  the  other  hand,  could 
occupy  a  large  extent  both  in  Doppler  and  in  range.  Figure  1  shows  the  signal  received  from  an 
HFSWR  in  the  rainge-Doppler  domain.  In  this  figure,  the  radar  look  direction  was  at  boresight. 
The  radar  transmits  a  frequency  modulated  interrupted  continuous  waveform  (FMICW)  with  a 
sweep-bandwidth  of  125  kHz  at  a  radar  frequency  of  5.8  MHz.  Two  clusters  of  ionospheric 
interference  can  be  seen:  (i)  from  around  100  km  to  over  150  km  (associated  with  reflections 
from  the  E-layer),  and  (ii)  from  about  210  km  to  beyond  300  km  (associated  with  reflections 
from  the  F-layer).  The  Doppler  extent  of  the  ionospheric  interference  was  a  few  Hertz. 
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Figure  1.  Typical  ionospheric  qlutter  in  range-Doppler  (velocity)  space. 


An  important  characteristic  of  ionospheric  clutter  is  that  it  is  range  specific,  and  its  statistical 
properties  may  be  estimated  more  reliably  using  time-averaging  of  the  sample  from  the  same 
range  bin. 

Co-channel  communications  interference  is  the  dominant  source  of  interference  at 
nighttime.  Occasionally  it  could  be  present  during  day-light  hours.  Figure  2  shows  the  received 
signal  from  one  of  the  receiver  channels  of  the  Cape  Bonavista  HFSWR  at  about  1800  hour  in 
the  Doppler-range  domain  after  which  the  transmitter  was  turned  off.  There  are  six  plots  in 
Figure  2,  corresponding  to  six  contiguous  CITs  of  20-second  duration.  It  can  be  seen  that  during 
the  last  two  CITs  the  noise  level  had  risen  significantly.  In  this  figure,  the  Doppler  axis  has  200 
points,  and  the  Doppler  increment  was  0.049  Hz.  Hence  the  Doppler  domain  was  about  9.8  Hz, 
centred  about  0  Hz.  Figure  3  shows  the  bearing-Doppler  profile  of  the  signal  received  by  the 
radar  for  a  particular  range  bin.  The  plots  in  Figure  3  were  obtained  by  forming  beams  pointing 
at  different  directions  using  the  antenna  array  samples.  These  plots  show  that  there  was  a  distinct 
directional  characteristic  to  this  interference.  More  importantly,  it  is  observed  that  the  directional 
characteristic  of  the  interference  was  more  or  less  range  independent.  This  permits  the  use  of 
range-averaging  to  extract  the  statistical  parameters  of  the  interference  components  required  for 
cancellation. 

2.  METHODOLOGY 

The  objectives  of  this  study  are  to  evaluate  quantitatively  the  accuracy  of  the  bearing 
estimate  in  the  presence  of  interference  and  to  determine  whether  there  is  a  difference  in  the 
accuracy  of  the  bearing  estimate  resulting  from  changing  the  order  of  operations  among  Doppler 
processing,  interference  cancellation  and  digital  beamforming.  On  the  basis  of  the  above  analysis 
we  shall  recommend  an  optimal  order  of  processing. 

To  achieve  the  above  objectives,  the  following  steps  will  be  followed:  (i)  select  the 
appropriate  data  set,  (ii)  inject  synthesized  target  samples  in  the  experimental  data,  (iii)  select  the 
appropriate  interference  cancellation  scheme,  (iv)  carry  out  Doppler  processing,  interference 
cancellation  and  digital  beamforming  in  various  orders  of  execution,  and  (iv)  compare  the 
accuracy  of  the  target  bearing  estimate  and  the  improvement  achieved  by  various  processing 
schemes. 

2.1  Interference  cancellation  algorithms 

Interference  suppression  or  cancellation  algorithms  exploit  the  directional  characteristics 
of  the  interfering  signals  by  applying  a  set  of  optimum  weights  to  the  radar  receiver  samples  such 
that  the  contribution  of  the  interference  components  to  the  array  output  is  minimized.  These 
techniques  were  originally  developed  for  communications  systems  in  which  both  the  desired 
signal  and  the  interference  are  present  over  an  extended  period  of  time.  For  radars,  the  desired 
signal  is  the  echo  from  a  target,  which  is  present  only  for  a  short  duration  (the  length  of  the 
transmit  pulse).  The  interference,  on  the  other  hand,  would  be  present  over  a  much  longer  period 
of  time. 


There  are  many  different  interference  suppression  algorithms  available.  In  this  study  we 
shall  consider  a  particular  method  based  on  coherent  side-lobe  cancellation  (CSLC)  [3,4]  because 
this  method  is  most  suited  for  the  signal  processing  hardware  architecture  implemented  in  the 
radar  developed  for  surveillance  of  the  east  coast  of  Canada. 

The  general  schematic  diagram  of  the  coherent  sidelobe  cancellation  algorithms  is  shown 
in  Figure  4.  There  are  N  main-array  antennas,  arranged  in  a  linear  array,  which  are  used  to  form 
beams  pointing  at  various  directions.  In  addition  there  are  M  auxiliary-array  antennas,  which  may 
or  may  not  be  collinear  with  the  main-array  antennas,  used  exclusively  for  deriving  the 
parameters  used  in  the  interference  cancellation  schemes.  Each  antenna  element  is  coupled  to  an 
independent  coherent  receiver  from  which  time  series  for  each  range  bin  can  be  obtained.  The 
interference  cancellation  is  carried  out  for  each  range  bin. 

The  derivation  of  the  basic  cancellation  algorithm  is  well  documented.  We  shall  quote 
references  and  describe  the  algorithm  below. 

Referring  to  Figure  4,  the  signals  from  the  main  and  auxiliary  antennas  for  a  given  range 
bin  are  represented  by  vectors  x  =  [x0,  x„  ...  xN.,]T  and  x  =  [y0,  y„  ...,  yM_,]T,  respectively,  where 
the  subscript  is  the  channel  (or  antenna  element)  index,  and  the  superscript,  T,  denotes  transpose. 
We  shall  refer  to  x  and  x  as  the  main-array  snapshot  and  auxiliary-array  snapshot,  respectively. 

An  estimate  of  the  interference  component  in  x  is  given  by  z,  where 


and 


—  [Zo>  Zp  •  .  •  >  Zn-j]1 


(1) 


z 


n 


(2) 


Vector  w,,  is  the  set  of  optimal  weights  for  the  n-th  main-array  element.  The  superscript  H 
denotes  complex  conjugate  transpose.  Vector  Xc  is  the  auxiliary  array  snapshot  corresponding  to 
the  main  array  snapshot  from  which  the  interference  estimate  is  to  be  subtracted.  It  is  not  to  be 
confused  with  x  which  are  used  to  compute  the  covariance  matrix  and  the  correlation  vectors. 

The  interference  estimates  are  subtracted  from  the  main-array  snapshot,  and  the  resulting 
signal  vector  is  x  =  [x0,  x„  .  .  .,  xN.,]T,  where 


-  W  H  V 
n  c 


(3) 


The  modified  snapshot  is  then  fed  to  the  beamformer,  yielding: 
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Figure  4.  Schematic  diagram  of  an  interference  cancellation  algorithm  . 


u(0)  =  vHx 


where  v  is  the  main  array  beam  steering  vector  given  by: 


v  =  [1,  e 


.2nd  .  Q 

i - sin0 

J  X 


■  2ZE(NIl)dsine 

•  ,  e  *  ]T. 


Here,  X  is  the  radar  wavelength,  d  is  the  antenna  element  spacing  and  0  is  the  beam  steering 
angle. 
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In  (5)  the  following  conventions  are  adopted:  (a)  element  n=0  is  the  left-most  antenna  and  (b) 
0  is  measured  clockwise  from  boresight. 

The  optimum  weight  vector  w„  is  obtained  from  the  product  between  the  inverse  of  the 
covariance  matrix  Ryy  and  the  correlation  vector,  r„: 

W  =  R,:1  r  («) 

n  yy  n 


The  covariance  matrix  and  the  correlation  vectors  are  defined  as: 


R„  -  Efc  XH]  (7) 


and 


rn  =  m  <] 


(8) 


for  n  =  0,1,2,..,N-1, 

respectively,  where  E  denotes  expected  value. 

The  correlation  vector  represents  the  correlation  between  the  elements  of  the  main  array  and 
those  of  the  auxiliary  array. 

Depending  on  the  type  of  interference  to  be  cancelled,  there  is  a  variety  of  ways  in  which 
the  covariance  matrix  and  the  correlation  vectors  can  be  computed.  The  most  common  ones  are 
(a)  time-averaging  and  (b)  range  averaging.  The  array  snapshots  are  samples  of  the  radar  returns 
corresponding  to  a  particular  range  bin.  Consequently  there  is  one  snapshot  for  each  transmitted 
waveform  per  range  bin.  Denoting  the  time  and  range  dependencies  of  the  snapshots  by  indices 
i  and  j,  respectively,  we  have: 

x(ij)  =  [Xfl(ij)>  x,(ij),  xN_j(ij)]T  (9) 


and 


y(U)  =  [  y0(iJ)>  y,(ij)>  •  •  •  • ,  yM-i(U)]T  (10) 
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2.1.1  Time-averaging. 


The  covariance  matrix  and  the  correlation  vectors  can  be  obtained  by  fixing  the  range 
index  j  and  averaging  over  a  number  of  preceding  snapshots,  i.e. 


y0(iJ)yo(iJ)  y0ftj)yroj) 


yoftjtyM-itfj) 


yi(M)yoGJ) 


yi(id)yM-i(U) 


yM-i(M)yo6J) 


••  ••  y^mu-iU 


(ii) 


and 


rn  =  [yodJXCio)  yiftjKCU)  ••  ••  yM-i 

^  i=l 


(12) 


where  I  is  the  number  of  time  samples  used  in  the  averaging. 


2.1.2  Range-averaging. 

Alternatively,  if  there  is  a  high  degree  of  correlation  in  the  interference  components 
among  snapshots  that  are  far  apart  in  range,  the  covariance  matrix  and  the  correlation  vectors  can 
be  obtained  by  fixing  the  snapshot  index  i  and  averaging  over  a  number  of  range  samples,  i.e. 


yotfjtyotfj)  y0(iJ)yi’(iJ) 


yiOd'tyoW) 


yiCd^M-iOj) 


yM_1(iJ)yo(ij) 


yM.iCy)yM-i(M) 


(13) 
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and: 


E 


n 


7E  MjK'Co)  yi(ijX(ij)  ••  ••  yM-i^oXcU) 

j  H 


T 


(14) 


where  J  is  the  number  of  range  samples  used  in  the  averaging. 

2.2  Order  of  execution. 

The  order  in  which  the  three  operations,  (a)  interference  suppression,  (b)  Doppler 
processing  and  (c)  digital  beam-forming,  are  executed  will  have  some  effect  on  the  over-all 
performance  in  terms  of  target  detection  and  bearing  estimate.  The  relative  merits  of  various 
orders  of  execution  will  also  influence  the  choice  of  methods  for  computing  the  covariance 
matrix  and  the  correlation  vectors.  We  shall  consider  each  case  and  determine  whether  they 
should  be  investigated  in  detail.  The  possible  combinations  are  the  following: 

(i)  Interference  suppression;  Doppler  processing;  digital  beamforming, 

(ii)  Interference  suppression;  digital  beamforming;  Doppler  processing, 

(iii)  Doppler  processing;  interference  suppression;  digital  beamforming, 

(iv)  Doppler  processing;  digital  beamforming;  interference  suppression, 

(v)  Digital  beamforming;  interference  suppression;  Doppler  processing, 

(vi)  Digital  beamforming;  Doppler  processing;  interference  suppression. 

Obviously,  cases  (i)  and  (ii)  should  yield  identical  results,  because  both  Doppler 
processing  and  digital  beamforming  are  performed  using  FFTs.  Since  the  FFT  is  a  linear 
operation,  the  order  of  execution  is  interchangeable.  The  same  argument  applies  to  cases  (iv)  and 
(vi).  Generally,  array  processing  will  alter  the  relative  amplitude  and  phase  relationships  among 
the  elements  and  may  produce  different  results  if  the  order  of  execution  is  interchanged  between 
interference  suppression  and  digital  beamforming.  This  would  be  the  case  if  the  main  array 
elements  are  also  used  in  the  computation  of  the  covariance  matrix.  For  the  interference 
cancellation  schemes  employed  in  this  study,  case  (ii)  and  case  (v)  should  yield  identical  results. 
The  reasons  follow. 

In  case  (v)  the  main-array  snapshot  is  first  multiplied  by  the  beam-steering  vector  v 
to  form  a  beam  in  direction: 


u'(6)  =  vHx 


(15) 


The  result  can  then  be  considered  as  the  output  of  a  single  antenna  element,  and  interference 
cancellation  as  defined  in  (1)  and  (2)  can  be  applied,  yielding: 
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(16) 


Z'  =  WH  Xc  =  (RyyV)\ 


and 


tL  =  X  (uV  =  y(vHx)*  (17) 

Subtracting  the  interference  estimate  z'  from  (14)  yields 

n'  -  z'  =  vH  (x-z)  =  u  (18) 


In  subsequent  discussions,  we  shall  refer  to  the  class  of  interference  suppression  schemes  that  are 
based  on  coherent  sidelobe  cancellation  as  "interference  cancellation".  For  this  class,  there  are 
only  two  distinct  cases  that  need  to  be  considered: 

Interference  suppression  in  the  time-domain  (i.e.  before  Doppler  processing), 
followed  by  Doppler  processing  and  beamforming; 

Interference  suppression  in  the  Doppler  domain  (i.e.  Doppler  processing  followed 
by  interference  cancellation  and  digital  beamforming). 

For  each  of  the  above  cases,  there  are  the  options  of  using  time-averaging  or  range¬ 
averaging  to  compute  the  covariance  matrix  and  the  correlation  vectors.  For  the  case  of 
interference  cancellation  in  Doppler  domain,  there  is  a  concern  in  the  appropriateness  of  using 
time-averaging  to  calculate  the  covariance  matrix  and  correlation  vectors  because  a  large  number 
of  pulses  is  required  to  form  the  Doppler  spectrum.  Hence,  a  significant  amount  of  time  would 
have  elapsed  if  a  number  of  spectra  is  to  be  used,  making  decorrelation  of  the  samples  highly 
likely. 

We  shall  investigate  the  following  cases: 

(a)  perform  interference  cancellation  in  the  time  domain  (i.e.  pulse  by  pulse)  using  time¬ 
averaging  for  the  computation  of  the  correlation  estimates. 

(b)  perform  interference  cancellation  in  the  time  domain  using  range-averaging  for  the 
computation  of  the  correlation  estimates. 

(c)  perform  interference  cancellation  in  the  Doppler  domain  (i.e.  Doppler  by  Doppler)  using 
the  range-averaging  method  to  compute  the  correlation  estimates. 
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2.2.1  Various  configurations 


Figure  5  shows  a  schematic  diagram  of  the  basic  coherent  sidelobe  cancellation  algorithm. 
For  a  particular  implementation  of  the  interference  cancellation,  the  appropriate  samples  from 
each  antenna  element  (both  main  and  auxiliary)  are  placed  in  a  data  buffer.  A  segment  of  the  data 
in  the  buffer,  starting  from  sample  #1!  and  ending  at  sample  #I2,  is  used  for  the  calculation  of  the 
covariance  matrix  and  correlation  vectors.  The  inverted  covariance  matrix  and  the  correlation 
vectors  are  used  to  derive  the  weight  vector  w„.  To  perform  cancellation  on  a  particular  data 
sample  in  the  buffer,  say  the  i-th  sample  of  the  n-th  main  antenna  element  (xn(i)),  the 
corresponding  sample  of  auxiliary  array  data  vector  is  multiplied  with  the  weight  vector  w„  and 
subtracted  from  the  main  antenna  sample  xn(i). 

The  calculation  of  the  covariance  matrix  and  the  correlation  vectors  is  essentially  a 
process  of  extracting  the  correlation  properties  of  the  interference  process  among  antenna 
elements.  Hence  we  shall  refer  to  these,  collectively,  as  correlation  estimates.  There  is  a  variety 
of  ways  to  obtain  the  correlation  estimates,  depending  on  what  type  of  data  are  placed  in  the  data 
buffers  in  Figure  5  and  how  they  are  utilized. 

For  interference  processes  that  are  isolated  in  range,  such  as  ionospheric  clutter,  the  use 
of  time  averaging  of  samples  from  a  fixed  range  bin  to  get  the  correlation  estimates  is  more 
appropriate. 

The  disadvantage  in  employing  time  averaging  for  the  correlation  estimates  is  that  it 
requires  a  number  of  snapshots  to  obtain  the  results.  If  the  time  between  the  estimation  process 
and  the  application  of  the  results  exceeds  the  decorrelation  time  of  the  interference  process,  poor 
cancellation  will  result.  Consequently,  it  is  necessary  to  employ  a  sliding  window  so  that  the  time 
lag  between  the  estimate  and  cancellation  process  is  minimized. 

Another  disadvantage  is  that,  since  data  samples  used  for  the  correlation  estimates  also 
come  from  the  range  bin  to  which  cancellation  is  to  be  applied,  there  is  a  possibility  that  the 
estimates  includes  the  target  component,  thereby  cancelling  some  of  the  target  component  as  well. 
This  cannot  be  avoided.  However,  in  most  cases  the  interference  energy  is  significantly  higher 
than  the  target  energy  and  improvement  can  still  be  obtained. 

The  following  is  a  list  of  different  configurations  whose  interference  suppression 
performance  will  be  investigated. 

(a)  Time-domain  cancellation  with  sliding-time-window  averaging 

Figure  6  shows  the  type  of  data  and  boundaries  of  the  data  used  for  the  calculation  of  the 
correlation  estimates  for  the  case  of  time-domain  cancellation  with  sliding-time-window 
averaging.  In  this  case,  the  data  to  be  stored  in  the  data  buffer  are  time  samples  (either  individual 
channel  or  beam-formed)  of  a  fixed  range  bin  over  a  number  (NFFT  +  NAVG)  of  consecutive 
pulses.  It  is  assumed  that  a  Doppler  processor  follows  operating  on  NFFT  cancelled  samples.  For 
each  time  sample  to  be  cancelled,  NAVG  preceding  samples  are  used  for  the  computation  of  the 


12 


data  buffer 


ixc 


13 


Figure  5.  Schematic  diagram  of  a  coherent  side-lobe  canceller. 
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Figure  6.  Time-domain  cancellation  with  sliding-time-  window  averaging. 
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Figure  7.  Time-domain  cancellation  with  fixed-range-window  averaging. 
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Figure  8.  Time-domain  cancellation  with  sliding-range-window  averaging. 
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Figure  9.  Doppler-domain  cancellation  with  fixed-range-window  averaging. 
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Figure  10.  Doppler-domain  cancellation  with  sliding-range-  window  averaging. 
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correlation  estimates  (hence  the  number  NFFT  +  NAVG). 

(b)  Time-domain  cancellation  with  fixed-range-window  averaging. 

Figure  7  shows  the  type  of  data  and  the  boundaries  of  the  data  used  for  calculating  the 
correlation  estimates  for  the  case  of  "Time-domain  cancellation  with  fixed  window  range 

averaging".  In  this  case,  the  data  to  be  stored  in  the  data  buffer  are  time  samples  (either 

individual  channel  or  beam-formed)  of  a  given  snapshot  for  all  the  range  bins.  The  last  NAVG 
range  samples  are  used  to  calculate  the  correlation  estimates. 

(c)  Time-domain  cancellation  with  sliding-range-window  averaging. 

Figure  8  shows  the  type  of  data  and  the  boundaries  of  the  data  used  for  calculating  the 
correlation  estimates  for  the  case  of  "Time-domain  cancellation  with  sliding-range-window 

averaging".  In  this  case,  the  data  to  be  stored  in  the  data  buffer  are  time  samples  (either 

individual  channel  or  beam-formed)  of  a  given  snapshot  for  all  the  range  bins.  The  data  samples 
used  for  the  calculation  of  the  correlation  estimates  are  divided  into  two  parts,  separated  by  a 
guard  zone  on  either  side  of  the  sample  to  be  cancelled.  The  guard  zones  are  used  to  ensure  that 
the  target  component  will  be  excluded  from  the  correlation  estimates. 

(d)  Doppler-domain  cancellation  with  fixed-range-window  averaging. 

A  similar  scheme  as  described  in  (b)  may  be  used  to  perform  cancellation  on  Doppler 
samples.  Figure  9  shows  the  type  of  data  and  the  boundaries  of  the  data  used  for  calculating  the 
correlation  estimates  for  the  case  of  "Doppler-domain  cancellation  with  fixed-range-window 
averaging".  In  this  case,  the  data  to  be  stored  in  the  data  buffer  are  samples  of  a  particular 
Doppler  of  a  given  snapshot  for  all  range  bins.  A  number  (NAVG)  of  samples  (usually  near  the 
maximum  range)  are  used  to  calculate  the  correlation  estimates.  The  results  are  used  to  cancel 
the  interference  component  of  the  samples  at  other  range  bins. 

(e)  Doppler-domain  cancellation  with  sliding-range-window  averaging. 

A  similar  scheme  as  that  described  in  (c)  may  be  used  to  perform  cancellation  on  Doppler 
samples.  Figure  10  show  the  type  of  data  and  the  boundaries  of  the  data  used  for  calculating  the 
correlation  estimates  for  the  case  of  "Doppler-domain  cancellation  with  sliding-range-window 
averaging". 

2.3  Types  of  interference 

In  this  study  we  shall  consider  the  cancellation  of  ionospheric  clutter  and  co-channel 
communications  interference  only. 

2.3.1  Ionospheric  clutter. 

Ionospheric  clutter  is  present  in  groups  of  range  bins  corresponding  to  altitudes  of  the  E 
and  F  layers  of  the  ionosphere.  One  would  not  expect  that  there  is  much  correlation  between  the 
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snapshots  from  two  range  bins  that  are  far  apart.  Consequently,  it  may  not  be  appropriate  to  use 
range-averaging  to  calculate  the  covariance  matrix  and  the  correlation  vectors.  On  the  other  hand, 
if  the  range  bins  used  in  the  calculation  of  the  covariance  matrix  and  the  correlation  vectors  are 
reasonably  close  to  the  range  bin  for  which  interference  cancellation  is  being  applied,  there  may 
be  sufficient  correlation  among  these  samples  to  permit  cancellation.  The  analysis  will  be 
confined  to  those  range  bins  which  contain  ionospheric  clutter,  and  the  range  bins  used  for 
computing  the  covariance  matrix  and  the  correlation  vectors  are  those  close  to  the  range  bin  for 
which  cancellation  is  to  be  applied. 

2.3.2  Co-channel  communications  interference. 

For  communications  interference,  the  snapshots  contain  the  samples  of  the  interference 
at  time  instants  that  are  separated  by  only  a  small  time  interval.  For  examples,  if  the  waveform 
repetition  frequency  (WRF)  of  the  radar  is  100  Hz,  then  the  snapshot  of  the  maximum  range  is 
separated  from  that  of  the  minimum  range  by  only  10  milli-seconds.  Since  the  interference  signal 
is  usually  comprised  of  continuous  signals  that  last  over  a  time  much  longer  than  the  radar’s 
waveform  repetition  interval  (WRI),  it  would  be  present  in  all  range  bins,  and  there  would  be  a 
high  degree  of  correlation  among  the  samples  taken  from  range  bins  that  are  far  apart. 
Consequently  the  correlation  properties  can  be  estimated  using  samples  from  one  group  of  range 
bins  and  the  resulting  weights  can  be  used  to  cancel  the  interference  at  other  ranges. 

2.4  Performance  criteria 

Some  care  must  be  exercised  in  defining  the  performance  criteria.  Our  objectives  are  to 
determine  (a)  whether,  and  to  what  degree,  the  various  interference  suppression  schemes 
introduce  errors  in  the  target-bearing  estimate,  and  (b)  which  of  the  interference  suppression 
schemes  is  optimum.  The  evaluation  must  be  carried  out  under  a  common  set  of  conditions  such 
as  SINR. 

2.4.1  Accuracy  of  the  bearing  estimate. 

The  measure  of  bearing-estimate  accuracy  is  relatively  straightforward.  We  use  the  mean 
bearing  error  and  its  standard  deviation  as  measures.  Synthesized  target  samples  with  given 
Doppler  and  angle-of-arrival  can  be  generated  and  embedded  in  experimental  data  that  contain 
interference.  The  data  can  then  be  processed  with  the  various  interference  suppression  schemes. 
Since  the  true  target  bearing  is  known,  the  error  in  target  bearing  can  be  evaluated.  The  mean 
error  and  standard  deviation  are  obtained  over  a  number  of  trials  with  independent  sets  of  data 
samples. 

There  is  one  problem,  and  that  is  the  interference  level  is  not  constant,  either  in  time,  or 
in  Doppler.  This  means  that  a  SINR  cannot  be  specified  independent  of  the  data.  Figure  1 1  shows 
the  Doppler  spectrum  of  the  radar  returns  from  a  particular  range  bin  (range  =  60  km).  The  data 
were  taken  from  file  AIR_ONE  (3  November  1995),  and  the  CIT  interval  was  approximately  28 
seconds.  Interference  is  observed  in  a  band  of  Doppler  frequencies  centred  at  about  -0.7  Hz; 


16 


Doppler  Frequency  (Hz) 

Figure  11.  Spectrum  of  radar  returns  with  a  high  level  of  communications  interference. 


however,  there  is  a  significant  amount  of  fluctuation  of  the  interference  level  over  the  Doppler 
band.  Ultimately  the  SINR  that  counts  is  the  one  specified  at  the  target  Doppler  shift.  Obviously, 
the  SINR  would  depend  on  the  target  Doppler. 

The  approach  we  adopted  is  the  following.  The  data  to  be  analyzed  were  divided  into 
blocks,  each  of  which  comprises  a  number  of  time  samples  in  a  CIT  interval.  The  angle-of-arrival 
of  the  interference  component  was  estimated  from  the  data.  Synthesized  target  samples  were 
generated  at  a  number  of  angle-of-arrivals  and  Doppler  shifts.  The  bearing  estimates  were 
compared  among  various  cases  that  include  (a)  no  interference  cancellation  and  (b)  with 
interference  cancellation.  The  analysis  was  repeated  for  a  number  of  SINR  values  (e.g.  +5, 0,  and 
-5  dB,  etc).  This  implies  that  the  actual  target  magnitude  is  dependent  on  the  interference  level 
at  the  specified  Doppler  in  a  CIT.  The  mean  and  standard  deviation  of  the  error  in  bearing 
estimated  were  obtained  over  the  number  of  CITs  analyzed. 

2.4.2  Improvement  factor 

To  arrive  at  a  figure-of-merit  for  interference  suppression  schemes  requires  some  thought. 
One  possible  performance  criteria  is  the  improvement  factor  usually  defined  as  the  ratio  between 
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the  SINR  after  and  before  the  application  of  interference  suppression.  This  requires  that  the  signal 
magnitude  be  known.  In  defining  the  signal  level,  it  is  preferable  that  a  fixed  reference  be  used. 
However,  this  is  not  always  possible  since  it  is  data  dependent. 

As  stated  earlier,  the  composite  signal  of  an  HFSWR  comprises  sea-clutter,  external  noise, 
ionospheric  clutter  and  possibly  co-channel  interference.  At  close-in  ranges,  the  sea  clutter  is 
dominant.  However,  the  Bragg  lines  could  still  be  a  significant  part  of  the  total  radar  return  even 
at  very  long  range.  Ionospheric  clutter  is  present  in  specific  ranges  and  has  essentially  a  low-pass 
spectrum.  The  same  problem  discussed  in  Section  2.4.1  is  also  present  here. 

In  radar,  because  of  the  constant-false-alarm-rate  (CFAR)  requirement,  the  detection 
performance  depends  not  only  on  the  SINR  at  the  target  Doppler  and  range,  but  also  on  the 
relative  magnitudes  of  the  noise  and  interference  in  neighbouring  range  and  Doppler  bins.  A  high 
level  of  noise  and  interference  in  neighbouring  Doppler  and  range  bins  results  in  a  high  detection 
threshold  setting,  thereby  lowering  the  probability  of  detection. 

There  are  two  problems  in  using  the  conventional  definition  of  improvement  factor.  These 
are  (a)  the  SINR  is  dependent  on  the  target’s  Doppler  shift  and  (b)  the  improvement  depends  not 
only  on  how  much  the  interference  is  suppressed,  but  also  on  how  much  the  signal  is  reduced 
by  the  interference-suppression  algorithm.  Since  there  is  no  a  priori  knowledge  regarding  the 
target  Doppler  and  angle  of  arrival,  an  input  SINR  can  only  be  defined  with  respect  to  a 
particular  set  of  data. 

For  the  interference  suppression  schemes  considered  in  this  study,  the  amount  of  reduction 
in  the  target-signal  magnitude  is  inversely  proportional  to  the  angular  distance  between  the  target 
and  the  interference.  The  reason  is  that  the  interference-suppression  algorithm  attempts  to  place 
a  null  on  the  antenna  pattern,  in  the  direction  of  the  interference.  If  the  bearing  of  the  target 
coincides  with  that  of  the  interference,  it  will  be  completely  cancelled,  thereby  yielding  no 
improvement  or  even  degradation.  On  the  other  hand,  if  the  bearing  of  the  target  is  well  separated 
from  that  of  the  interference  (say,  more  than  one  beamwidth),  there  is  generally  a  negligible 
amount  of  reduction  in  the  target-signal  magnitude.  Consequently,  a  reasonable  alternative 
definition  of  the  improvement  factor  is  as  follows: 


Average  interference  power  before  suppression 
Average  interference  power  after  suppression 


(19) 


where  the  average  is  over  Doppler,  range  or  both. 

The  above  definition  of  the  improvement  factor  is  only  valid  for  the  case  where  the  target 
bearing  is  sufficiently  separated  from  that  of  the  interference.  In  evaluating  the  improvement 
factor,  some  times  it  is  necessary  to  exclude  other  components  such  as  sea  clutter  because  these 
may  be  the  dominant  component  in  the  signal  on  which  the  interference  suppression  algorithm 
may  not  have  any  effect.  More  will  be  said  about  this  in  Section  3  which  describes  the  data 
analysis. 
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2.5 


Test  data 


Data  collected  in  experimental  trials  conducted  at  Cape  Bonavista  and  Cape  Race  have 
been  analyzed  to  determine  their  suitability  for  this  study. 

2.5.1  Cape  Bonavista  data 

Most  of  the  data  obtained  from  Cape  Bonavista  trials  were  collected  during  daylight 
hours.  The  few  data  sets  that  were  collected  at  nighttime  contain  system  artifacts  that  make  them 
unsuitable  for  this  study.  In  particular,  data  from  four  out  of  the  eight  channels  contain  an 
injected  calibration  tone.  A  considerable  amount  of  effort  would  be  required  to  remove  this 
component  from  the  data  which  is  necessary  for  the  case  of  time-domain  interference 
cancellation.  Consequently,  the  Cape  Bonavista  data  will  not  be  used  in  this  study. 

2.5.2  Cape  Race  data 

Three  sets  of  data  collected  from  Cape  Race  will  be  used  in  this  study.  Table  1  tabulates 
the  essential  parameters  of  these  data  files. 


File 

Name 

Date 

Time 

Frequency 

Waveform 

PRF 

(Hz) 

No.  of 
channels 

AIR_ONE 

3-May 

1995 

00:39 

5.81 

FMICW 

9.01 

10 

AIRJTVE 

27-Oct 

1994 

11:57 

5.81 

FMICW 

9.18 

10 

AIR_SEVEN 

28-Oct 

1994 

18:06 

5.81 

FMICW 

9.18 

10 

Table  1.  Cape  Race  data  sets  used  in  this  study. 


The  Cape  Race  data  has  the  following  characteristics.  First,  the  ten  receiver  channels  were 
connected  to  ten  sub-arrays  comprising  four  receive-antenna  elements.  The  element  separation 
is  22.22  m  which  is  equalled  to  half  a  wavelength  at  6.75  MHz.  Hence  the  effective  spacing  of 
the  sub-arrays  is  88.88  m.  This  results  in  grating  lobes  in  the  receive-antenna  pattern.  However, 
because  the  performance  of  the  various  processing  schemes  is  evaluated  on  a  relative  basis,  it  is 
not  considered  to  be  a  serious  problem.  Second,  there  exist  range-ambiguity  responses  arising 
from  the  sub-optimal  pulse-compression  algorithm  used  in  the  Cape  Race  radar.  These  anomalies 
could  corrupt  the  covariance  matrix  and  correlation  vector  calculations.  Fortunately  the  range- 
ambiguity  responses  were  isolated  in  range  and  their  locations  were  known.  Consequently  those 
samples  corresponding  to  the  range-ambiguity  response  were  excluded  in  the  calculations.  In  the 
analysis  that  follows,  Channels  1  and  10  were  used  as  auxiliary  antennas.  Channels  2  to  9  forms 
an  eight-element  main  array. 
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Data  file  AIR_ONE  was  collected  around  midnight.  These  data  had  a  high  external  noise 
level  as  well  as  communications  interference.  This  file  was  used  to  comparatively  evaluate  the 
various  interference  suppression  schemes  and  their  effects  on  the  accuracy  of  bearing  estimate. 

Date  file  AIR_FIVE  was  collected  at  mid  day.  The  data  in  a  number  of  range  bins  contain 
strong  ionospheric  clutter.  This  file  was  used  to  test  algorithms  for  possible  ionospheric  clutter 
suppression. 

Data  file  AIR_SEVEN  was  collected  at  dusk.  The  data  in  this  file  contained  echoes  of  a 
controlled  target  as  well  as  some  communications  interference.  This  file  was  used  to  compare 
various  interference  suppression  schemes  under  practical  conditions. 

3.  ANALYSIS  RESULTS. 

The  Cape  Race  Radar  had  ten  receiver  channels.  Unless  otherwise  stated,  the  results 
presented  in  this  section  were  obtained  under  the  following  conditions. 

(i)  Main  array  comprised  channels  2-9. 

(ii)  Auxiliary  array  comprised  channels  1  and  10 

(iii)  No.  of  snapshots  used  for  the  evaluation  of  correlation  properties  =  64. 

3.1  Communications  interference. 

3.1.1  Case  where  the  communications  interference  is  dominant. 

Data  file  AIR_ONE  listed  in  Table  1  were  collected  in  a  trial  on  May  1995  around  mid¬ 
night.  At  this  time,  the  external  noise  level  was  very  high.  Figure  12  shows  the  range-Doppler 
profile  of  the  radar  samples  in  AIR_ONE  at  boresight  over  one  CIT.  It  shows  that  sea  clutter 
(Bragg  lines)  was  visible  out  to  only  about  70  km.  The  external  noise  floor  was  at  approximately 
-90  dB,  and  there  was  interference  around  the  centre  of  the  Doppler  band  that  was  more  or  less 
range  independent. 

Figure  13  shows  the  results  of  applying  time-domain  cancellation  to  the  same  data  as  in 
Figure  12,  using  sliding-time-window  averaging  to  calculate  the  correlation  estimates.  Each  time 
sample  was  cancelled  using  the  correlation  estimates  computed  from  64  preceding  snapshots.  It 
shows  that  the  results  were  worse  than  in  Figure  12  where  no  cancellation  was  applied. 

Figure  14  shows  the  results  of  applying  time-domain  cancellation  to  the  data,  using  fixed- 
range-window  averaging  to  calculate  the  correlation  estimates.  The  time  sample  in  each  range 
bin  was  cancelled  using  the  correlation  estimates  computed  from  64  range  samples  positioned 
near  the  maximum  range.  The  interference  level  had  been  reduced. 
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Figure  15  shows  the  results  of  applying  time-domain  cancellation  to  the  data,  using 
sliding-range-window  averaging  to  calculate  the  correlation  estimates.  Sixty  four  range  samples 
near  the  sample  to  be  cancelled  were  used  to  calculate  the  correlation  estimates.  The  guard  zones 
contain  one  range  bin.  This  was  sufficient  for  the  Cape  Race  Radar  because  its  range  resolution 
was  one  range  samples.  The  results  show  that  the  interference  level  had  been  suppressed. 

Figure  16  shows  the  results  of  applying  Doppler-domain  cancellation  to  the  data,  using 
fixed-range-window  averaging  to  calculate  the  correlation  estimates.  The  correlation  estimates 
were  obtained  the  same  way  as  that  for  Figure  14,  except  that  the  time  samples  were  first 
Doppler  processed.  Cancellation  was  applied  to  the  Doppler  samples,  and  the  interference  level 
was  reduced. 

Figure  17  shows  the  results  of  applying  Doppler-domain  cancellation  to  the  data,  using 
sliding-range-window  averaging  to  calculate  the  correlation  estimates.  The  correlation  estimates 
were  obtained  the  same  way  as  that  for  Figure  15,  except  that  Doppler  samples  were  used. 
Reduced  interference  level  was  obtained. 

The  relative  merits  of  the  various  processing  schemes  can  be  evaluated  by  comparing  the 
averaged  interference  levels  after  cancellation  was  applied.  Because  the  communications 
interference  was  dominant  and  was  more  or  less  range  independent,  the  Doppler  spectra  can  be 
averaged  over  all  the  range  bins.  The  results  are  shown  in  Figure  18a  through  18f,  for  the  spectra 
in  Figures  12-17,  respectively.  It  can  be  seen  in  Figure  18a  that,  before  cancellation  was  applied, 
the  interference  occupied  a  Doppler  domain  from  -2  Hz  to  about  +1  Hz. 

The  result  in  Figure  18b  shows  that  the  performance  was  worse  using  time-domain 
cancellation  with  sliding-time-window  averaging.  The  reasons  for  this  performance  degradation 
will  be  discussed  in  Section  4.2. 

Results  in  Figures  18c  through  18f  showed  that,  qualitatively,  time-domain  cancellation 
with  sliding-range-window  averaging  yields  the  best  performance.  The  improvement  factor  may 
be  derived  from  the  results  of  Figures  18a.  The  two  peaks  around  zero  Doppler  were  those  of  the 
ground  clutter  at  close  ranges  and  the  Bragg  components.  The  spectra  in  Figure  18  can  be 
averaged  over  the  Doppler  domain,  with  the  samples  around  zero  Doppler  excluded.  The  results 
are  shown  in  Table  2. 

3.1.2  Case  where  the  communications  interference  is  not  dominant. 

Data  file  AIR_SEVEN  listed  in  Table  1  was  collected  in  a  trial  on  28  October,  1994  at 
dusk.  At  that  time,  the  external  noise  level  had  begun  to  increase,  and  there  was  some 
communications  interference.  Figure  19  shows  the  range-Doppler  profile  of  the  radar  samples 
over  one  CIT.  It  shows  that  the  noise  floor  was  at  approximately  -1 15  dB.  The  centre  portion  of 
the  Doppler  spectrum  was  occupied  by  the  sea  clutter  which  extended  out  to  beyond  300  km.  In 
the  Doppler  region  between  2.3  Hz  and  4.0  Hz,  there  was  a  band  of  interference  which  appears 
in  all  range  bins.  The  level  of  this  interference  was  substantially  below  that  of  the  sea  clutter. 
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Figure  17.  Range-Doppler  profile  of  radar  returns  with  Doppler-domain  cancellation,  using  sliding-range-window  averaging. 


Magnitude  (dB)  Magnitude  (dB)  Magnitude  (dB) 


FILE:  AIRONE;  FREQ.  =  5.81  MHz;  WRF  =  9.01  Hz;  CIT  #  30;  CIT  =  28  sec. 
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Figure  18.  Averaged  spectra  before  and  after  interference  cancellation. 
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Avg. 

Interference 

level 

Improvement 

Factor 

(AIR.ONE) 

(a) 

No  interference  cancellation 

-88.41 

(b) 

time-domain  cancellation  using  sliding-time- 
window  averaging 

-79.64 

-8.77 

(c) 

time-domain  cancellation  using  fixed-range-window 
averaging 

-92.39 

3.98 

(d) 

time-domain  cancellation  using  sliding-range- 
window  averaging 

-96.07 

7.66 

(e) 

Doppler-domain  cancellation  using  fixed-range- 
window  averaging 

mm 

2.04 

(f) 

Doppler-domain  cancellation  using  sliding-range- 
window  averaging 

-94.41 

6.00 

Table  2:  Comparison  of  Improvement  Factor  among  various  schemes,  for  the  case  where 

the  communications  interference  is  dominant. 


At  the  time  this  data  set  was  collected,  there  was  a  Challenger  601  aircraft  flying  towards 
the  Cape  Race  radar  along  boresight.  The  speed  of  the  aircraft  was  about  178  knots.  The 
corresponding  target  Doppler  fell  in  the  region  occupied  by  the  interference,  and  the  aircraft  was 
not  detected.  Figures  20  through  24  show  the  range-Doppler  profiles  of  the  radar  samples  after 
interference  cancellation  was  applied. 

Figures  20,  21  and  22  are  for  cases  of  time-domain  cancellation  with  (i)  sliding-time- 
window  averaging,  (ii)  fixed-range-window  averaging  and  (iii)  sliding-range-window  averaging, 
respectively.  These  results  show  that  time-domain  cancellation  failed  to  suppress  the 
communications  interference. 

Figures  23  and  24  are  for  cases  of  Doppler-domain  cancellation,  with  (i)  fixed-range- 
window  averaging  and  (ii)  sliding-range-window  averaging,  respectively.  These  results  show  that 
Doppler-domain  cancellation  is  capable  of  suppressing  the  communications  interference,  and  the 
aircraft  target  emerged  after  cancellation. 

The  spectra  in  Figures  20  through  24  were  averaged  over  all  range  bins,  and  the  results 
are  shown  in  Figure  25.  In  cases  where  time-domain  cancellation  was  used,  the  result  was 
generally  worse  than  the  case  where  no  cancellation  was  used.  In  the  case  of  time-domain 
cancellation,  with  sliding-range-window  averaging,  artifacts  in  the  response  can  result  (see  Figure 
25d). 
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Figure  20.  Range-Doppler  profile  of  radar  returns  with  time-domain  cancellation  using  sliding-time-window  averaging. 


Figure  21.  Range-Doppler  profile  of  radar  returns  with  time-domain  cancellation  using  fixed-range- window  averaging. 
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Figure  24.  Range-Doppler  profile  of  radar  returns  with  Doppler-domain  cancellation,  using  sliding-range-window  averaging. 


Magnitude  (dB)  Magnitude  (dB)  Magnitude  (dB) 


To  get  an  estimate  of  the  improvement  factor,  the  spectra  in  Figure  25  were  averaged  over 
the  Doppler  domain.  In  this  case,  we  need  to  exclude  the  sea  clutter  because  it  is  not  part  of  the 
interference  against  the  target.  Hence  the  spectra  were  averaged  over  the  Doppler  bins  from  1  Hz 
to  4.59  Hz  only.  The  results  are  shown  in  Table  3. 

For  the  case  where  the  communications  interference  was  not  dominant,  Doppler-domain 
cancellation  with  sliding-range-window  averaging  yielded  the  best  result. 


Avg. 

Interference 

level 

Improvement 

Factor 

(AIR_SEVEN) 

(a) 

No  interference  cancellation 

-113.95 

- 

(b) 

time-domain  cancellation  using  sliding-time- 
window  averaging 

-80.06 

-33.89 

(c) 

time-domain  cancellation  using  fixed-range- 
window  averaging 

-99.14 

-14.81 

(d) 

time-domain  cancellation  using  sliding-range- 
window  averaging 

-111.32 

-2.63 

(e) 

Doppler-domain  cancellation  using  fixed-range- 
window  averaging 

-119.31 

5.36 

(f) 

Doppler-domain  cancellation  using  sliding-range- 
window  averaging 

-120.06 

6.11 

Table  3:  Comparison  of  Improvement  Factor  among  various  schemes,  for  the  case  where 

the  communications  interference  is  not  dominant. 


3.1.3  Effect  of  interference  cancellation  on  accuracy  of  bearing  estimate. 

The  analysis  presented  in  this  subsection  is  intended  to  give  an  indication  of  the  relative 
accuracy  of  the  bearing  estimate  between  two  different  interference  cancellation  procedures.  It 
is  not  intended  as  an  absolute  quantitative  measure  of  accuracy.  The  difficulty  in  obtaining  an 
absolute  quantitative  measure  of  accuracy  arises  from  the  fact  that  it  depends  on  many  factors 
which  are  not  known  a  priori.  These  include  the  spectral  and  directional  characteristics  of  the 
interference  signals,  and  the  separation  of  the  target,  both  in  Doppler  and  in  bearing,  from  the 
interference.  Nevertheless,  we  shall  make  an  attempt  to  get  some  answers  with  respect  to  a  set 
of  specific  conditions. 
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Three  methods  of  processing  radar  data  were  studied: 


(a)  No  Sidelobe  Cancellation  (NO_SLC),  i.e.  Doppler  Processing  only, 

(b)  Doppler  Processing  followed  by  Sidelobe  Cancellation  (DP_SLC),  and 

(c)  Sidelobe  Cancellation  followed  by  Doppler  Processing  (SLC_DP). 

The  emphasis  was  on  target  detectability  and  accuracy  of  direction  estimates. 

The  processing  parameters  and  definitions  are  briefly  described  in  what  follows. 

(i)  Each  Doppler  spectrum  has  I,  where  I  =  256,  amplitudes.  It  was  obtained  as  the  inverse 
FFT  of  the  output  of  an  array  element  for  a  specified  range  bin  in  I  consecutive  data 
blocks. 

(ii)  The  frequency  bins  in  the  Doppler  spectra  are  identified  with  i  =  -1/2  to  i  =  1/2,  and  the 
Doppler  frequency  at  the  i-th  bin  is  given  by: 

f  =  PRF  (20) 

I 

where  PRF  is  the  pulse  repetition  frequency.  An  approaching  target  produces  a  positive  Doppler 

frequencies. 

(iii)  For  the  radar  frequency  of  5.81  MHz,  the  Bragg  lines  are  at  frequencies  of  -0.2459  and 
0.2459  Hz  (Doppler  bins  i  =  -7  and  i  =  7),  respectively. 

(iv)  Each  directional  spectrum  has  K,  where  K  =  256,  amplitudes.  It  was  obtained  as  the 
inverse  FFT  of  a  K-point  sequence  obtained  by  augmenting  the  N  amplitudes  of  a 
processed  main  array  snapshot  with  K-N  zeros.  The  processing  included  Doppler 
processing  and,  if  needed  adaptive  sidelobe  cancellation. 

(v)  The  bins  in  the  K-point  FFT  direction  spectra  are  labelled  with  k  =  -K/2  to  k  =  K/2.  The 
direction  at  the  k-th  bin  is  given  by: 


0k  =  —  (beamwidth)  (21) 

K 

or 

0k  =  sin'1/^-)  (radians)  (22) 

Kd 

where  N  is  the  number  of  main-array  elements,  d  is  the  element  spacing  and  X  is  the  radar 
wavelength.  All  directions  are  measured  clockwise  from  boresight. 
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by: 


A  direction  in  beamwidths,  0bw,  is  related  to  the  corresponding  direction  in  radians,  0rad, 


0 


bw 


rad  * 


(23) 


(vi)  A  target  direction  estimate,  denoted  by  0,  is  an  estimated  obtained  from  a  K-point  FFT 
directional  spectrum.  It  is  the  direction  of  the  spectral  peak  closest  to  the  true  target 
direction  0.  The  error  in  the  estimate  is  defined  as: 


A0  =  0  -  0 


(24) 


A  smaller  A0  corresponds  to  a  more  accurate  direction  estimate. 

(vii)  Let  x’  be  a  target-absent  array  snapshot  in  a  specified  Doppler  bin.  The  interference-plus- 
noise  power  in  this  bin  is  calculated  as: 


P  +  P  = 

rl  rN 


N 


(25) 


Let  x"  be  the  target-only  array  snapshot  at  the  same  Doppler  bin.  The  signal  power  in  the 
bin  is  calculated  as: 


The  composite  snapshot,  with  both  target  and  interference  present  is  defined  as: 


x  =  x'  + 


x" 


(26) 


(27) 


the  input  SINR  in  x  is  calculated  as: 


SINRj 


?s 

Pi  +Pn 


(28) 
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(viii)  Output  SINR  in  a  processed  snapshot  x  is  calculated  as: 


SINR0 


(29) 


where  Xj  and  Xj.  are  the  components  of  x  parallel  and  orthogonal  to  the  array  steering  vector 
v(§),  respectively.  A  larger  SINR0  corresponds  to  a  higher  detection  probability. 

Given  a  number  of  pairs  of  {2L?»2L}’  the  output  SINR  is  calculated  as: 


SINR0  = 


Ave(  |X|  |2) 
Ave(  |x^  |2) 


(30) 


Average  values  were  used,  because  they  smoothed  out  the  effect  of  interaction  between  the 
injected  and  interference  signals. 

A  preliminary  examination  of  the  data  showed  that  the  night-time  interference  was  present 
in  all  range  and  Doppler  bins.  For  this  data  set,  the  interference  peaked  at  Doppler  bin  i  =  34, 
corresponding  to  1.2189  Hz.  The  direction  of  the  peak  interference  was  observed  to  vary  with 
time  between  direction  bins  k  =  32  and  k  =  100.  After  Doppler  filtering,  the  highest  peak  was 
close  to  the  direction  bin  k  =  66  in  all  Doppler  bins.  Bins  k  =  32,  100  and  66  corresponded  to 
1.0,  3.12  and  2.06  beamwidths,  respectively,  measured  clockwise  from  the  boresight  direction. 

The  above  characteristic  regarding  peak  position  is  consistent  with  the  fact  that  the 
ionosphere  is  highly  dynamic.  This  directional  characteristic  is  essentially  independent  of  Doppler 
and  range. 

File  AIR_ONE  was  modified  by  the  injection  of  a  target  with  a  specified  range,  Doppler 
frequency  and  direction.  One  hundred  and  five  non-overlapping  sets,  each  containing  I,  where 
I  =  256,  consecutive  array  snapshots  with  the  specified  range  were  extracted  from  the  modified 
file.  From  each  set,  an  array  snapshot  x  with  the  specified  Doppler  frequency  was  generated, 
using  the  methods  described  in  Section  3.13a,  b  and  c.  This  xwas  then  used  to  calculate  a  target 
direction  estimate  Sand  a  vector  pair  {>»,&},  where  and  4  were  the  components  of  x parallel 
and  orthogonal  to  the  array  steering  vector  v(^,  respectively.  The  105  values  of  b,  and  Xj_  were 
then  used  to  calculate  an  average  error  in  direction  estimate,  Ave  { A()} ,  a  standard  deviation  of 
direction  estimates,  Std{§},  and  an  output  SINR,  SINR0.  To  facilitate  the  comparison  of  results, 
all  directions  are  expressed  in  beamwidths  (bw)  measured  clockwise  from  boresight. 

Table  4  contains  the  values  of  Ave{A0},  Std{0}  and  SINRq  calculated  with  sixteen  target 
directions  spaced  0.5  beamwidth  apart.  The  injected  target  signal  was  in  range  bin  50  and 
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Table  4:  Dependence  of  Ave(A9)  ,  Std(0)  and  SINR„  on  beam  direction  for  SINRj  =  5  dB.  The 

Doppler  bin  i  =  34  and  the  interference-plus-noise  power  Pi+PN=0.1975E-06. 
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Table  5:  Dependence  of  Ave(A0)  ,  Std(0)  and  SINR0  on  beam  direction  for  SINR(  =  0  dB.  The 

Doppler  bin  i  =  34  and  the  interference-plus-noise  power  Pj+PN=0.1975E-06. 
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Table  6:  Dependence  of  Ave(A0)  ,  Std(0)  and  SINR0  on  beam  direction  for  SINRj  =  -5  dB.  The 

Doppler  bin  i  =  34  and  the  interference-plus-noise  power  Pj+PN=0.1975E-06. 


Doppler  frequency  bin  i  =  34  and  its  SINR,  was  5.0  dB.  In  the  last  two  rows  are  the  averages 
and  standard  deviations  of  the  values  in  the  columns.  There  are  several  notable  observations. 
Firstly,  the  values  vary  with  directions  and  the  variation  is  largest  in  the  range  1.0<6<3.0  bw. 
Secondly,  the  direction  estimates  were  the  most  accurate  in  the  results  obtained  with  NO_SLC, 
less  accurate  with  SLC_DP  and  least  accurate  with  DP_SLC.  Thirdly,  the  average  SINRq  is  6.5 
dB  in  the  NO_SLC  results,  about  5.8  dB  in  the  SLC_DP  results,  and  less  than  3.0  dB  in  the 
DP_SLC  results.  Tables  6  and  7  contain  the  results  obtained  by  lowering  the  SINR,  to  0.0  and  - 
5.0  dB,  respectively.  From  the  average  and  standard  deviations  in  the  last  two  rows,  it  was 
observed  that  the  SLC_DP  method  produced  the  best  results. 

Table  7  contains  the  results  calculated  with  sixteen  Doppler  frequencies  spaced  16  bins 
apart.  The  injected  target  signal  was  in  range  bin  50  and  direction  0  =  0.  Its  SINR,  was  0.0  dB. 
With  the  exception  of  slightly  larger  values  of  Ave{A?)},  it  was  observed  that  the  SLC_DP 
produced  the  best  results  for  all  frequencies. 

Table  8  contains  the  results  calculated  with  four  different  range  bins  and  target  directions 
0  =  0.0,  0.31,  and  1.25  bw.  The  injected  target  signal  was  in  Doppler  bin  i  =  34  and  the  SINR, 
was  0.0  dB.  With  the  exception  of  larger  values  of  Ave{A§}  for  0  =  1.25,  it  was  observed  that 
the  SLC_DP  procedure  produced  the  best  results  for  all  range  bins. 

The  following  conclusions  are  based  on  the  results  in  Tables  4  to  8: 

(i)  The  NO_SLC  results  are  the  best  if  SINR,  =  5  dB  or  higher;  and 

(ii)  The  SLC_DP  results  are  the  best  if  SINR,  =  0  dB  or  lower  and  the  target  direction  0  is 

not  close  to  the  peak  direction  of  the  interference.  When  the  two  directions  are  close  the 

magnitude  of  Ave(A0}  is  larger  than  those  in  the  NO_SLC  and  DP_SLC  results. 

A  few  remarks  are  in  order.  Firstly,  the  SINR0  cannot  fall  below  a  certain  minimum  value; 
otherwise,  the  target  cannot  be  detected.  The  minimum  value  is  system  dependent  and  will  not 
be  discussed  here.  Based  on  this  consideration,  the  SLC_DP  method  provides  the  best  detection 
performance  when  the  SINR,  is  low. 

Secondly,  adaptive  sidelobe  cancellation  is  unnecessary  if  SINR,  is  high.  In  this  case,  one 
expects  the  results  to  be  better  than  those  obtained  with  SLC_DP  and  DP_SLC,  because  sidelobe 
cancellation  modifies  the  gain  patterns  of  the  main  array.  The  modification  changes  the  relative 
magnitudes  of  the  signals  and  noise  and  usually  adds  a  bias  to  the  target  direction  estimate  as 
well. 


Thirdly,  the  SLC_DP  results  are  usually  better  than  the  DP_SLC  results,  because  night¬ 
time  interference  is  produced  by  ionospheric  reflection  of  unwanted  target  signals.  The  ionosphere 
is  highly  dynamic  and  produces  rapid  fluctuations  in  arrival  directions  of  interference  signals.  The 
DP_SLC  method  averages  these  fluctuations  over  time  and  is  less  effective  than  the  SLC_DP 
method  which  works  on  the  instantaneous  interference  directions.  The  variation  of  the  parameters 
in  Table  4  with  direction  agrees  with  the  earlier  observation  that  the  interference  direction  is 
centred  near  2.06  beamwidth  off  boresight. 
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Table  7:  Dependence  of  Ave(A0)  ,  Std(0)  and  SINR0  on  Doppler  frequency  for  a  target  in  range  bin  50  and  bearing 

at  boresight  (SINR,  =  0  dB). 
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To  gain  some  insight  into  the  accuracy  of  the  bearing  estimate  of  a  real  target  after 
interference  cancellation,  the  data  segments  that  contain  the  echoes  of  the  Challenger  601  and 
communications  interference  were  analyzed.  These  data  segments  correspond  to  CIT  No.  70  to 
77  of  data  file  AIR_SEVEN.  In  these  CIT  intervals,  the  target  was  not  detected  because  of  the 
interference.  Figure  26a  shows  the  spectrum  of  the  data  in  CIT  No.7 1 .  The  Doppler  component 
corresponding  to  the  Challenger  is  submerged  in  interference.  Figure  26b  shows  the  spectrum  of 
the  same  data  after  Doppler-domain  cancellation  with  sliding-range-window  averaging  is  applied. 
The  interference  around  the  target  Doppler  has  been  suppressed  and  the  target  has  emerged. 

The  Doppler  shift  of  the  target  in  Figure  26  is  estimated  to  be  3.55  Hz  which  corresponds 
to  a  radial  velocity  of  178.2  knots.  The  angular  response  of  this  Doppler  component  is  obtained 
by  performing  digital  beam  forming  on  this  Doppler  component,  using  a  256-point  FFT.  The 
results  are  shown  in  Figure  27.  In  both  cases  (with  and  without  interference  cancellation),  the 
peak  response  of  the  target  is  at  0.39°.  This  agrees  well  with  the  ground-truth  data  of  the  target 
which  indicated  that  the  aircraft  flew  along  the  aperture  boresight  towards  the  radar. 

3.2  Ionospheric  clutter. 

There  are  more  constraints  on  the  cancellation  of  ionospheric  clutter  than  that  of 
communications  interference.  Ionospheric  clutter  is  a  self-generated  interference  and  not  present 
when  the  transmit  signal  or  a  reflective  medium  is  not  present.  Consequently,  it  is  not  appropriate 
to  use  range  averaging  to  calculate  the  correlation  estimates.  This  leaves  time-domain  cancellation 
with  sliding-time-window  averaging  as  the  only  viable  solution.  However,  as  demonstrated  in 
Section  3.1,  the  effectiveness  of  time-domain  cancellation  depends  on  how  quickly  the  correlation 
estimates  are  obtained  and  applied.  For  low  WRF  systems  such  as  the  Cape  Race  radar,  the  time 
between  snapshots  is  too  long  for  the  time-domain  cancellation  to  be  effective  since  it  takes  at 
least  several  tens  of  snapshots  to  obtain  reliable  correlation  estimates.  By  the  time  the  estimation 
process  is  completed,  several  seconds  would  have  elapsed,  rendering  the  cancellation  ineffective. 
Nevertheless,  if  the  WRF  is  sufficiently  high  such  that  the  time  between  snapshots  is  short,  time- 
domain  cancellation  will  yield  good  ionospheric  clutter  suppression. 

Since  the  WRF  of  the  radar  was  rather  low,  these  data  are  not  suitable  for  analysis  of 
time-domain  cancellation  schemes  for  ionospheric  clutter,  using  either  time  or  range  averaging. 
There  is  one  other  scheme  which  we  have  not  tried  up  to  this  point  that  might  work  when  the 
WRF  is  low.  That  is  "Doppler-domain  cancellation  with  Doppler  averaging".  The  rationale  is  as 
follows.  Assume  that  the  ionospheric  clutter  is  the  dominant  interference  in  the  returns  from  a 
number  of  range  bins  and  that  it  has  a  Doppler  spectrum  that  encompasses  a  moderate  number 
of  Doppler  bins.  If  the  clutter  is  coming  from  a  small  number  (smaller  than  the  number  of 
auxiliary  antennas)  of  bearings,  then  the  various  Doppler  components  associated  with  the 
ionospheric  clutter  will  have  similar  directional  characteristics.  And  the  covariance  matrix  and 
the  correlation  vectors  may  be  estimated  by  averaging  Doppler  samples.  To  examine  this 
approach  we  use  data  from  the  Cape  Race  radar. 
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(a)  Doppler  spectrum  with  no  cancellation. 


DOPPLER  FREQUENCY  (Hz) 


(b)  Doppler  spectrum  after  Doppler-domain  cancellation. 

Figure  26.  Spectra  of  data  containing  echoes  from  a  Challenger  601 . 
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Figure  27.  Angular  response  of  the  target  Doppler  component. 


Figure  28  shows  the  Doppler  spectrum  of  the  returns  at  boresight  from  range  bin  no.  200 
(about  240  km  from  Cape  Race).  The  data  were  taken  from  file  AIR_FIVE,  and  the  spectrum  was 
computed  using  a  2048-point  FFT.  This  FFT  size  corresponds  to  a  CIT  of  223  seconds  which  is 
compatible  with  ship  detection.  The  spectral  components  from  -0.4  Hz  to  about  +0.9  Hz  were 
those  of  the  ionospheric  clutter.  Since  the  WRF  of  the  data  was  9.18  Hz,  there  were  about  290 
Doppler  bins  in  this  region. 

A  simple  scheme  that  is  easy  to  implement  is  to  use  all  (except  the  one  to  which 
cancellation  is  to  be  applied  and  a  number  of  neighbouring  Doppler  bins  which  serves  as  a  guard 
zone)  the  Doppler  samples  in  the  region  (-0.4  Hz  to  +0.9  Hz)  to  compute  the  covariance  matrix 
and  the  correlation  vectors.  An  estimate  of  the  interference  derived  from  these  correlation 
estimates  is  subtracted  from  each  Doppler  sample  in  the  region.  Since  both  estimation  and 
cancellation  take  place  within  the  same  range  bin  at  the  same  time  instant,  the  problem  of  the 
long  elapsed  time  between  estimation  and  cancellation  is  circumvented. 

It  was  conjectured  that,  if  there  is  no  significant  variation  in  the  directional  characteristics 
of  the  ionospheric  clutter  within  a  CIT,  this  scheme  could  work.  Figure  29  shows  the  resulting 
spectrum  after  cancellation  was  applied  to  all  Doppler  samples,  using  the  scheme  described 
above.  The  ionospheric  clutter  was  suppressed  in  some  of  the  Doppler  components.  The  average 
improvement  factor  in  the  Doppler  region  (-0.4  Hz,  +0.9  Hz)  was  about  2  dB. 
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To  investigate  how  this  interference  cancellation  scheme  affects  a  target,  synthesized 
target  samples  were  generated  and  imbedded  in  the  time  series  before  Doppler  processing.  The 
target  bearing  was  at  boresight,  and  the  target  magnitude  was  -96  dB.  From  Figure  29,  it  can  be 
seen  that  the  suppression  was  greater  on  the  positive  side  of  the  Doppler  spectrum.  Hence  we 
chose  a  target  Doppler  of  0.15  Hz. 

Figure  30  shows  the  spectrum  of  the  time  series  with  the  imbedded  target  without 
interference  cancellation.  It  is  difficult  to  differentiate  the  target  from  the  interference.  Figure  31 
shows  the  resulting  spectrum  after  interference  cancellation  had  been  applied.  The  target  has  now 
emerged  from  the  clutter. 

From  the  above  analysis,  it  is  evident  that  the  coherent  sidelobe  cancellation  scheme  as 
implemented  for  co-channel  communications  interference  is  not  very  effective  against  ionospheric 
clutter.  The  reasons  are  as  follows. 

(i)  Ionospheric  clutter  could  come  from  several  directions,  and  these  directions  may  vary 
with  time.  Consequently,  the  auxiliary  antenna  array  must  have  a  sufficient  number  of 
elements  to  provide  the  required  degrees  of  freedom  to  cancel  the  spatially-distributed 
clutter. 

(ii)  In  the  ionospheric  clutter  region,  there  are  frequent  occurrences  of  meteoric  clutter.  The 
meteoric  clutter  alters  the  temporal  and  spatial  correlation  characteristics,  rendering  their 
estimation  inaccurate. 

(iii)  Ionospheric  clutter  is  isolated  in  range.  This  means  that  it  is  not  appropriate  the  employ 
range-averaging  to  estimate  the  correlation  properties  needed  for  coherent  sidelobe 
cancellation. 

(iv)  Ionospheric  clutter  is  time  varying  both  in  range  and  Doppler.  This  means  that  employing 
time  averaging  or  Doppler  averaging  to  estimate  the  correlation  properties  is  likely  to  be 
ineffective. 

It  may  be  possible  to  improve  the  performance  of  the  coherent  sidelobe  cancellation 
scheme  by  a  judicious  use  time-averaging,  range  averaging  and  Doppler  averaging  to  obtain  the 
correlation  properties.  Alternatively,  techniques  other  than  the  coherent  sidelobe  cancellation 
scheme  may  be  employed.  The  development  of  these  schemes  and  the  verification  of  their 
performance  is  beyond  the  scope  of  the  current  study.  Nevertheless,  it  is  important  to  gain  a  good 
understanding  of  the  behaviours  of  ionospheric  clutter  in  the  temporal,  spatial  and  spectral 
domains.  In  Appendix  A,  we  present  the  analysis  of  some  available  ionospheric  clutter  data  in 
terms  of  its  temporal,  spatial  and  spectral  behaviours. 

Based  on  the  results  presented  in  Appendix  A,  it  appears  that  the  ionospheric  clutter 
arising  from  overhead  reflection  off  the  E-region  is  less  of  a  problem  at  higher  frequencies  (about 
4  MHz)  where  ship  detection  is  concerned.  The  reason  for  this  is  that  the  magnitude  of  the 
ionospheric  clutter  from  the  E  region,  at  higher  frequencies,  is  significantly  lower  than  the  Bragg 
lines.  It  is  most  pronounced  during  the  morning  hours.  Since  the  apparent  range  of  the  E-region 
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Figure  30.  Doppler  spectrum  of  data  with  an  imbedded  synthetic  target 
with  no  cancellation. 
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Figure  31.  Doppler  Spectrum  of  data  with  an  imbedded  synthetic  target 
after  Doppler-domain  cancellation  using  Doppler  averaging. 


is  around  100  km,  vessels  with  a  moderately  large  RCS  would  still  yield  sufficient  signal-to- 
clutter  ratio  for  detection. 

At  longer  ranges  (250  km  to  350  km),  ionospheric  reflection  presents  a  bigger  problem. 
The  magnitude  of  the  F-region  ionospheric  clutter  at  higher  frequencies  is  higher  than  that  for 
the  E-region.  In  addition,  at  long  ranges,  the  target  magnitude  would  be  attenuated,  resulting  in 
a  much  lower  signal-to-clutter  ratio.  For  long-range  ship  detection  in  the  F-region  ionospheric 
clutter  zone,  the  radar  should  operate  at  frequencies  below  4  MHz. 

4.  CONCLUSIONS  AND  DISCUSSION. 

4.1  Conclusions 

From  analyses  presented  in  Section  3,  the  following  conclusions  can  be  drawn. 

(a)  For  interference  suppression  schemes  based  on  coherent  sidelobe  cancellation  techniques, 
there  are  only  two  distinct  cases  in  which  different  results  will  be  produced  as  a 
consequence  of  changing  the  order  of  execution  among  interference  suppression,  Doppler 
processing  and  digital  beam  forming.  These  are  (i)  interference  cancellation  before 
Doppler  processing  and  (ii)  interference  cancellation  after  Doppler  processing.  Changing 
the  order  of  execution  between  Doppler  processing  and  digital  beamforming  either  in  the 
time  domain  or  the  Doppler  domain  does  not  affect  the  outcome. 

(b)  For  the  case  where  co-channel  communications  interference  is  dominant,  the  optimum 
processing  procedures  are  as  follows. 

(i)  Interference  cancellation  on  each  time  sample  using  sliding-range-window 
averaging  for  computation  of  the  correlation  estimates; 

(ii)  Doppler  processing 

(iii)  Digital  beamforming. 

Steps  (ii)  and  (iii)  may  be  interchanged. 

(c)  For  the  case  where  co-channel  communications  interference  is  not  dominant,  the  optimum 
processing  procedures  are  as  follows. 

(i)  Doppler  processing; 

(ii)  Interference  cancellation  on  each  Doppler  sample  using  sliding-range-window 
averaging  for  calculation  of  the  correlation  estimates; 

(iii)  Digital  beamforming. 

Steps  (ii)  and  (iii)  may  be  interchanged. 

(d)  For  ionospheric-clutter  cancellation,  using  range  averaging  for  the  computation  of  the 
covariance  matrix  and  the  correlation  vectors  is  not  appropriate.  The  reason  for  this  is  that 
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ionospheric  clutter  is  present  in  a  small  number  of  range  bins  only,  and  there  is  very  little 
correlation  among  samples  from  range  bins  that  are  far  apart.  At  this  time,  the 
recommended  processing  procedures  is: 


(i)  Interference  cancellation  in  the  time-domain  (i.e.  before  Doppler  processing); 
using  sliding-time-window  averaging  for  the  computation  of  the  covariance  matrix 
and  the  correlation  vectors;  use  a  small  number  of  samples  for  the  average  to 
minimize  decorrelation; 

(ii)  Doppler  processing  and  beamforming  or  vice  versa. 

(e)  In  terms  of  the  effects  of  interference  suppression  on  the  accuracy  of  bearing  estimates, 

refer  to  Section  3.1.3. 

4.2  Discussion. 

The  effectiveness  of  any  interference  suppression  scheme  based  on  coherent  sidelobe 
cancellation  using  array  antennas  depends  on  several  factors.  The  first  is  that  there  must  be  some 
directivity  in  the  interference.  The  second  is  the  number  of  angles-of-arrival  of  the  interference 
signal.  This  dictates  the  number  of  degrees  of  freedom  the  system  must  have  to  yield  an  effective 
cancellation. 

The  third  is  that  the  time  between  the  estimation  of  the  correlation  properties  of  the 
interference  process  among  various  antenna  elements  and  the  application  of  cancellation  must  be 
much  shorter  than  the  decorrelation  time  of  the  interference  process.  For  example,  in  the  data 
analyzed  in  this  study,  the  WRF  was  about  9  Hz  and  thus  the  time  between  two  contiguous 
snapshots  is  about  1 10  msec.  If  the  interference  process  occupies  a  Doppler  domain  of  a  few  Hz, 
the  decorrelation  time  is  in  the  order  of  a  few  hundred  milli-seconds.  If  time-averaging  is  used 
to  calculate  the  correlation  estimates,  it  will  take  several  tens  of  snapshots  to  get  the  results.  By 
the  time  the  estimates  are  obtained,  several  seconds  would  have  elapsed,  and  poor  cancellation 
is  likely  to  result.  This  is  supported  by  the  analysis  results.  Table  9  compares  the  elapsed  time 
for  calculating  the  correlation  estimates  and  the  improvement  factor  (when  applied  to  data  set 
AIR_ONE)  among  three  different  schemes.  Progressively  improving  cancellation  results  are 
obtained  as  the  elapsed  time  decreases. 

Decorrelation,  in  itself,  does  not  present  problems  if  the  interference  originates  in  one 
spatial  location  only.  In  reality,  however,  interference  comes  from  many  directions.  Because  it 
takes  different  amount  of  time  for  the  various  interference  signals  to  arrive  at  the  antenna 
elements,  the  relative  amplitude  and  phase  relationships  among  various  antenna  elements  could 
change  significantly  with  time  due  to  decorrelation. 

Although  cancellation  in  the  Doppler-domain  does  not  always  yield  the  most 
improvement,  there  are,  however,  several  advantages  when  the  overall  performance  is  considered. 
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Correlation  estimation  scheme 

Elapsed 

time 

Improvement  factor 
(AIR  ONE) 

(dB) 

Sliding-time-window  averaging  (64  snapshots) 

6.9  sec. 

-8.99 

Fixed-range-window  averaging  (using  64  samples 
near  the  maximum  range) 

-110 

msec 

3.98 

Sliding-range-window  averaging  (64  range  samples) 

-hundreds 
of  psec. 

7.66 

Table  9:  Comparison  of  improvement  factor  as  a  function  of  elapsed  time  between 

correlation  estimates  and  cancellation. 


The  first  is  that  it  provides  a  means  for  isolating  the  various  interference  sources,  thereby 
allowing  the  algorithm  to  operate  on  weaker  interferences.  The  second  is  that  the  computation 
load  is  greatly  reduced.  For  example,  cancellation  is  not  required  for  Doppler  components  that 
are  outside  the  range  of  expected  target-Doppler  shifts.  In  addition,  the  available  time  for  the 
required  processing  is  measured  in  terms  of  tens  of  seconds  (depending  on  the  size  of  the  FFT), 
rather  than  milli-seconds,  as  in  the  case  of  pulse-by-pulse  cancellation.  The  third  is  that  it 
preserves  the  information  contained  in  the  data.  For  example,  time-domain  interference 
suppression  distorts  the  relative  magnitude  of  the  sea-clutter  components,  rendering  it  useless  for 
monitoring  the  ocean-surface  conditions.  Finally,  it  can  be  implemented  without  disturbing  the 
basic  signal-processing  structure.  Because  the  output  of  the  Doppler  processor  is  near  the  end  of 
the  signal  processing  chain,  the  interference-suppression  module  can  be  implemented  as  a  self- 
contained  module  inserted  between  the  Doppler  process  and  the  digital  beam  former. 

For  improved  low-altitude  aircraft  detection  against  communications  interference, 
suppression  in  the  Doppler  domain  is  the  most  appropriate  scheme.  The  reason  for  this  is  that 
the  range  of  interest  is  within  150  km  where  the  sea  clutter  is  dominant.  Doppler-domain 
suppression  permits  the  isolation  of  the  sea  clutter  from  the  interference. 

Concerning  the  suppression  of  ionospheric  clutter,  however,  there  is  insufficient 
information  at  present  to  draw  the  conclusion  as  to  whether  time-domain  or  Doppler-domain 
suppression  yields  a  better  performance.  Because  the  ionospheric  clutter  is  isolated  in  range,  it 
precludes  the  use  of  Doppler-domain  cancellation  using  range-averaging  for  the  correlation 
estimates.  From  the  limited  analysis  presented  in  Appendix  A,  it  appears  that  the  dominant 
ionospheric  clutter  comes  from  a  small  number  of  directions,  perhaps  3  to  4.  Consequently,  given 
a  sufficient  number  of  auxiliary  antenna  elements,  these  can  be  cancelled  effectively.  The 
difficulty  is  in  attempting  to  cancel  the  ionospheric  clutter  continuum.  The  clutter  can  originate 
from  a  large  number  of  directions  and  thus  requiring  a  prohibitively  large  number  of  auxiliary 
antennas.  Furthermore,  because  the  amplitude  of  the  clutter  continuum  is  significantly  lower  than 
the  specular  component,  it  would  be  difficult  to  estimate  its  correlation  properties  precisely  in 
the  presence  of  the  much  larger  specular  components. 


55 


There  is  one  other  property  that  we  should  examine.  That  is  the  true  angle-  of-arrival  of 
the  ionospheric  clutter  in  the  elevation  plane.  This  requires  a  planar  array  or  a  cross  array.  If  it 
can  be  established  that  the  angle-of-arrival  of  the  ionospheric  clutter  in  elevation  is  confined  to 
a  narrow  region,  then  it  may  be  possible  to  cancel  it  with  a  finite  number  of  auxiliary  antenna 
elements  configured  to  yield  nulls  in  elevation.  More  experiments  and  analyses  are  required 
before  a  definite  conclusion  can  be  drawn. 

5.  REFERENCES. 

[1]  R.  H.  Khan,  "Ocean-Clutter  Model  for  High-Frequency  Radar,"  IEEE  Journal  of  Ocean 
Engineering,  Vol.  62,  pp.  181-188,  April  1991. 

[2]  M.W.Y.  Poon,  R.H.  Khan  and  S.  Le-Ngoc,  "A  Singular  Value  Decomposition  based 
Method  for  Suppressing  Ocean  Clutter  in  High  Frequency  Radar",  IEEE  Trans  on  Signal 
Processing,  Vol.41,  No. 3,  March  1993,  pp.1421-1425. 

[3]  Howells,  "Intermediate  Frequency  Sidelobe  Canceller,"  U.S.  Patent  3  202  990,  August 
1965. 

[4]  S.P.  Applebaum,  "Adaptive  Arrays,"  IEEE  Transactions  on  Antennas  and  Propagation, 
AP-24,  September  1976,  pp.  585-598. 

[5]  L.  Petrie,  private  communication,  Petrie  Telecommunication,  October  1999. 

[6]  W.R.  Pigott  and  K.  Rawer,  Ed.,  "URSI  Handbook  of  Ionogram  Interpretation  and 
Reduction,"  Amsterdam,  New  York:  Elsevier,  1961. 

[7]  J.  Capon,  R.J.Greenfield  and  R.J.  Kolker,  "Multidimensional  Maximum  Likelihood 
processing  of  a  large  aperture  seismic  array",  Proceedings  of  the  IEEE,  Vol.55,  No.2, 
February  1967,  pp.  192-211. 


6.  ACKNOWLEDGEMENTS. 

The  authors  thank  Dr.  A.M.  Ponsford  of  Raytheon  Systems  Canada  Limited  for  useful 
discussions.  Special  thanks  are  due  Mr.  Len  Petrie  of  Petrie  Telecommunications  Limited  for 
providing  the  interpreted  ionogram  data.  This  work  was  supported  partly  by  Raytheon  Systems 
Canada  Limited  and  the  Canadian  Department  of  National  Defence. 


56 


7.  Appendix  A:  A  cursory  look  at  the  temporal,  spatial  and  spectral  behaviours  of 
ionospheric  clutter. 

To  develop  effective  algorithms  to  combat  ionospheric  clutter  in  HFSWR,  a  thorough 
understanding  of  its  behaviour  is  essential.  Ideally,  one  would  like  to  observe  ionospheric  clutter 
in  terms  of  its  strength,  Doppler  spread,  spread  in  range,  angles  of  arrival,  temporal  correlation, 
spatial  correlation.  In  addition,  it  is  also  desirable  to  observe  any  change  in  behaviour  as 
functions  of  time-of-day,  season  as  well  as  radar-system  parameters  such  as  frequency, 
bandwidth,  waveform,  etc.  This  would  entail  a  long  period  of  data  collection  and  analysis.  In  this 
appendix  we  present  a  preliminary  examination  of  the  characteristics  of  ionospheric  clutter.  This 
examination  is  limited  in  scope  because  of  the  limited  amount  of  available  data.  It  is  understood 
that  these  observations  were  site-specific  and  that  only  a  limited  number  of  frequencies  were 
used.  Consequently,  no  generalization  of  these  results  should  be  made.  Nevertheless,  some  insight 
may  be  gained  from  these  observations  that  could  tell  us  the  probable  effectiveness  of  a  particular 
algorithm.  For  example,  if  it  is  observed  that  the  ionospheric  clutter  originates  from  many 
directions  (both  in  azimuth  and  elevation),  then  sidelobe  cancellation  techniques  would  require 
too  many  degrees  of  freedom  to  be  effective. 

Data  collected  from  the  April  1995  iceberg-tracking  trial  at  Cape  Bonavista  were  used. 
The  Cape  Bonavista  radar  employed  a  simple  pulsed  waveform  with  the  following  parameters. 

Frequency  =4.1  MHz. 

PRF  =  12.5  Hz. 

Waveform  =  50  psec,  raised  cosine  shaping. 

Transmit  array  =  5  element  log-periodic  array,  Beamwidth  *  100°. 

Receive  array  =  8  element  doublet  array,  element  spacing  =  38.46  m. 

(the  two  monopole  antennas  in  a  doublet  was  spaced  nominally 
a  quarter  wavelength  apart  at  3.9  MHz) 

A.  1  Range-Doppler  views  of  ionospheric  clutter. 

Figure  A.  1  shows  the  range-Doppler  characteristics  of  the  ionospheric  clutter  observed  by 
the  Cape  Bonavista  HFSWR.  The  data  are  taken  from  the  April  1995  iceberg  tracking  trial.  The 
radar  parameters  are  shown  in  the  figure.  The  look  direction  was  at  boresight  (Beam  No.  1,  110° 
true).  Since  the  radar  frequency  was  4. 1  MHz,  the  Bragg  lines  were  at  fB  =  ±0.206  Hz,  where 
fB  is  the  Bragg  resonance  frequency  given  by: 


where  g  =  9.81  m/sec2  is  the  gravity  acceleration,  and  X  is  the  radar  wavelength. 
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Figure  A.  1  A  typical  range-Doppler  profile  of  ionospheric  clutter. 


The  Bragg  lines  are  shown  in  the  figure  as  two  relatively  large  Doppler  components  around  fB 
that  extend  out  to  beyond  300  km.  The  CIT  interval  was  327.68  seconds  (4096  pulses),  a  value 
suitable  for  ship  detection.  At  approximately  270  km,  a  spectral  components  can  be  observed 
between  the  two  Bragg  lines  with  a  range  extent  of  about  20  to  30  km.  This  component  is 
attributed  to  the  backscatter  from  the  ionosphere,  probably  the  F  region.  Depending  on  the  season 
and  the  time-of-day,  the  altitude  of  the  F  region  may  vary  between  250  km  to  beyond  400  km. 
Since  the  waveform  used  in  this  radar  had  a  nominal  bandwidth  of  20  kHz,  the  range  resolution 
cell  was  approximately  7.5  km.  Hence  the  thickness  of  the  region  that  gives  rise  to  this 
ionospheric  clutter  component  was  relatively  small,  perhaps  no  more  than  15  to  20  km.  For  the 
data  set,  ionospheric  reflection  around  100  km  range  (nominally  the  E  region)  was  observed  only 
intermittently,  and  its  magnitude  was  small  compared  with  those  observed  at  ranges  beyond  250 
km  (nominally  the  F  region). 

The  Doppler  spectrum  of  the  270-km  range  bin  in  Figure  A.l  is  shown  in  Figure  A.2.  The 
ionospheric  clutter  occupied  a  Doppler  region  between  ±0.4  Hz.  It  can  be  observed  that,  even 
though  there  was  a  peak  in  the  spectrum,  it  was  not  significantly  higher  than  its  neighbours. 
Figure  A.3  shows  the  spectrum  of  another  time  series,  taken  at  1400  hour,  of  a  range  bin  in 
which  the  peak  ionospheric  clutter  was  observed.  A  dominant  spectral  peak  was  observed  at 
0.067  Hz,  with  a  magnitude  close  to  20  dB  greater  than  the  next  highest  peak.  This  peak  could 
be  the  result  of  specular  reflection  from  some  specific  regions  of  the  ionosphere.  Thus  the  peak 
of  the  ionospheric  clutter  drifted  in  range  by  36  km  in  about  two  hours. 

Examining  additional  data,  it  was  observed  that  the  location  of  the  peak  ionospheric 
clutter  in  both  range  and  Doppler,  as  well  as  its  magnitude,  varied  as  a  function  of  the  time-of- 
day.  Figure  A.4  shows  the  apparent  height  (range)  of  the  spectral  peak  as  a  function  of 
Newfoundland  local  time.  It  shows  that  the  ionospheric  clutter  was  at  the  greatest  range  (apparent 
height)  in  the  afternoon  and  was  at  the  closest  range  just  before  midnight.  In  addition,  it  was 
observed  that,  during  daylight  hours,  there  was  a  greater  short-term  variation  (i.e.,  from  one  CIT 
to  the  next)  in  the  range  of  the  peak  ionospheric  clutter. 

The  magnitude  of  the  peak  ionospheric-clutter  component  also  varied  as  a  function  of  the 
time-of-day.  Figure  A. 5  shows  the  relative  magnitude  of  the  spectral  peak  as  a  function  of 
Newfoundland  local  time.  It  shows  that  the  ionospheric  clutter  was  at  the  lowest  level  around 
noon  and  the  highest  after  dusk  but  before  midnight.  Again,  during  daylight  hours,  a  greater 
short-term  variation  in  the  magnitude  of  the  peak  was  observed.  The  magnitudes  spanned  a  range 
of  about  30  dB  between  midday  and  midnight. 

No  significant  variation  in  the  Doppler  shift  of  the  peak  component  was  observed  as  a 
function  of  the  time-of-day.  Figure  A.6  shows  that  Doppler  shift  of  the  peak  ionospheric 
component  as  a  function  of  Newfoundland  local  time.  The  Doppler  shift  of  the  peak  varied 
between  ±0.2  Hz.  There  was,  however,  a  slight  bias  towards  the  positive  Doppler. 

The  assumption  that  the  observed  ionospheric  clutter  was  that  reflected  from  the  F-region 
was  corroborated  by  calculations  carried  out  by  Petrie  [5,6].  He  produced  prediction  of  the 
apparent  height  of  the  reflection  and  the  relative  signal  strength  as  a  function  of  time  based  on 
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DOPPLER  FREQUENCY  (Hz) 


RELATIVE  MAGNITUDE  OF  PEAK  ECHO  (dB)  APARENT  HEIGHT  OF  PEAK  ECHO  (km) 


NEWFOUNDLAND  LOCAL  TIME  (Hr) 


Figure  A.4  Variation  of  apparent  height  of  ionospheric  clutter. 


NEWFOUNDLAND  LOCAL  TIME  (Hr) 


Figure  A.5  Variation  of  peak  magnitude  of  ionospheric  clutter. 
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the  ionosonde  data  of  30  April  1995.  These  results  are  shown  in  Figures  A.7  and  A.8, 
respectively.  Good  qualitative  agreements  between  these  and  the  experimental  data  are  observed. 
In  both  the  prediction  and  the  radar  data,  the  maximum  and  minimum  height  occurred  at  around 
1 600  and  2000  hours,  respectively.  The  actual  maximum  range  observed  by  the  radar  was  greater 
by  a  few  tens  of  km;  however,  this  could  partly  be  attributed  to  the  normal  tolerance  of  the 
model  and  partly  to  the  fact  that  the  prediction  assumed  a  vertical  incidence  of  the  reflected 
signal.  The  direction  of  the  ionospheric  clutter  may  be  at  some  angle  off  the  vertical  axis,  thereby 
yielding  a  larger  slant  range. 

There  was  also  qualitative  agreement  between  the  predicted  and  observed  relative  signal 
strength  of  the  ionospheric  clutter.  Both  were  observed  to  be  at  a  minimum  at  1300  hour.  The 
range  of  variation  of  the  predicted  signal  strength  was  much  less  than  that  in  Figure  A.5; 
however,  in  Figure  A.5,  the  magnitude  was  that  of  the  peak  Doppler  component  and  not  the 
actual  signal  power  of  the  total  ionospheric  clutter. 

A.2  Temporal-spatial  behaviour  of  ionospheric  clutter. 

The  temporal  behaviour  of  ionospheric  clutter  may  be  observed  by  examining  the  time 
series  of  a  range  bin  that  contain  ionospheric  clutter.  Figure  A.9  shows  the  in-phase  waveform 
of  such  a  time  series  over  a  time  interval  of  20  minutes.  It  can  be  observed  that  a  significant 
level  of  ionospheric  clutter  was  present  in  about  half  of  the  time.  With  a  typical  CIT  interval  of 
200  seconds  for  ship  detection,  this  implies  that  ship  detection  will  be  affected  for  3  or  4  CITs. 
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Figure  A. 9  In-phase  waveform  of  a  time  series  containing  ionospheric  clutter. 


A  number  of  questions  are  of  interest  in  regard  to  the  development  of  effective 
ionospheric  clutter  suppression  schemes.  The  first  is  "what  are  the  angles  of  arrival  of  the 
ionospheric  clutter  and  how  do  they  vary  with  time?"  We  already  know  that  the  direction  of  the 
ionospheric  clutter  is  probably  near  the  vertical  direction.  What  is  not  known  is  the  exact 
elevation  angle  and  the  number  of  directions  in  azimuth.  To  determine  the  true  angles  of  arrival 
of  a  signal,  one  needs  a  two-dimensional  array.  Since  we  had  only  data  collected  from  a  linear 
array,  we  could  only  examine  the  apparent  angles-of-arrival  in  azimuth. 

To  determine  the  angles-of-arrival  of  the  ionospheric  clutter  in  azimuth,  we  took  the  range 
samples  from  all  the  array  elements  and  used  the  Capon  [7]  method  to  compute  the  angular 
response.  The  Capon  method  has  a  better  angular  resolution  than  the  FFT  approach  for  the  same 
number  of  antenna  elements.  For  the  data  in  Figure  A.9,  we  noticed  that  the  ionospheric  clutter 
began  to  be  visible  at  time  t  >  250  sec.  We  examined  256  contiguous  snapshots  starting  at  time 
t=255  sec.  Since  the  PRF  of  the  radar  was  12.5  Hz,  the  snapshots  were  80  msec  apart,  and  total 
time  interval  was  20.48  sec.  Figure  A.  10  shows  the  angular  responses  of  the  radar  for  a  particular 
range  (315  km)  as  a  function  of  snapshot  number.  It  can  be  observed  that,  during  this  20-second 
time  interval,  there  were  three  prominent  peaks,  located  at  nominal  azimuths  of  -68°,  5°  and  26°, 
respectively.  The  mean  and  standard  deviation  of  these  three  peaks,  over  the  256  snapshots,  were 
-67.83°,  +4.75°,  25.91°  and  1.72°,  1.95°,  and  1.91°,  respectively.  Hence  the  angles-of-arrival  of 
the  dominant  peaks  did  not  vary  appreciably  over  the  20-second  time  interval. 

The  two  peaks  near  +5°  and  +26°  were  not  unusual  because  they  could  still  come  from 
a  near-vertical  direction  but  slightly  off  boresight  in  azimuth.  The  peak  near  -68°,  however,  was 
puzzling  because  of  the  large  apparent  azimuthal  angle.  A  linear  array  cannot  provide  both 
azimuthal  and  elevation  information  on  a  signal  simultaneously.  If  the  signal  is  from  the  vertical 
direction,  then  the  apparent  azimuthal  angle  (as  observed  by  a  linear  array)  should  be  close  to 
zero.  However,  since  the  linear  array  employed  in  the  Cape  Bonavista  was  comprised  of 
monopole  doublets,  it  is  possible  that  increased  mutual  coupling  effects  could  distort  the 
elemental  pattern  at  large  squinting  angles.  Further  investigation  will  be  required  to  arrive  at  an 
explanation  for  this  phenomenon. 

Another  question  of  interest  is  whether  there  is  any  correlation  among  the  angles-of-arrival 
of  the  various  Doppler  components  that  constitute  the  ionospheric  clutter  spectrum.  This  is  of 
interest  because,  if  the  directional  characteristics  of  one  Doppler  component  are  similar  to  those 
for  another,  then  the  correlation  properties  estimated  from  one  Doppler  component  may  be  used 
to  cancel  another  Doppler.  Figure  A.l  1  shows  the  Doppler  spectra  of  a  time  series,  at  a  range  of 
323  km,  for  six  contiguous  CITs  of  163.8-second  duration.  Strong  ionospheric  clutter  components 
centred  about  zero  Doppler  were  evident  in  the  spectra.  In  addition,  several  spectral  peaks  were 
observed  whose  magnitudes  were  significantly  higher  than  the  clutter  continuum. 

Figure  A.  12  compares  the  angular  responses  of  the  largest  Doppler  components  in  the 
spectrum  of  Figure  A. lib.  In  this  figure,  the  strongest  Doppler  component  was  at  +0.134  Hz, 
followed  by  the  ones  at  -0.006  Hz  and  -0.195  Hz,  respectively.  We  compared  the  angular 
responses  of  these  three  Doppler  components,  together  with  one  (at  -0.0916  Hz)  selected  from 
the  clutter  continuum.  It  can  be  observed  that  the  correlation  of  the  angular  response  was  not  very 
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Figure  A.  10  Variation  of  angles-of-arrival  of  ionospheric  clutter  with  time. 
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Figure  A.  11  Spectra  of  ionospheric  clutter  in  consecutive  CITs. 
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CAPE  BONAVISTA 


Figure  A.  12  Angular  response  of  different  Doppler  components  of  the  ionospheric  clutter  in  a  particular  range  bin 


good  between  the  four  Doppler  components.  This  implies  that  there  is  no  guarantee  that 
correlation  characteristics  estimated  from  one  Doppler  component  can  be  used  to  cancel  another 
Doppler  component  effectively. 
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